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^ PBEFACE. 



The present publication is an essay that was sent in (De- 
cember, 1887) to compete for the Smith's Prizes at Cambridge. 
To the onlooker it is always a mournful thing to see what 
he considers splendid abilities or opportunities wasted for lack of 
knowledge of some paltry common-place truth. Such is in the 
main my feeling when considering the neglect of the study of 
Quaternions by that wonderful corporation the University of Cam- 
bridge. To the alumnus she is apt to appear as the leader in all 
branches of Mathematics. To the outsider she appears rather as 
the leader in Applied Mathematics and as a ready welcomer of 
other branches. 

If Quaternions were simply a branch of Pure Matheujatics we 
could understand why the study of them was almost confined 
to the University which gave birth to them, but as the truth is 
quite otherwise it is hard to shew good reason why they have not 
struck root also in Cambridge. The prophet on whom Hamilton's 

^ mantle has fallen is more than a mathematician and more than a 

natural philosopher — he is both, and it is to be noted also that he 

-. is a Cambridge man. He has preached in season and out of 

season (if that were possible) that Quaternions are especially useful 

* in Physical applications. Why then has his Alma Mater turned a 

~ deaf ear ? I cannot believe that she is in her dotaofe and has lost 

^ . . . . 

her hearing. The problem is beyond me and I give it up. 

But I wish to add my little efforts to Prof. Tait's powerful 
advocacy to bring about another state of affairs. Cambridge is the 
prepared ground on which if anywhere the study of the Physical 
applications of Quaternions ought to flourish. 
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VI PREFACE. 

When I sent in the essay I had a^ faint misgiving that per- 
chance there was not a single man in Cambridge who could 
understand it without much labour — and yet it is a straightfor- 
ward application of Hamilton's principles. I cannot say what 
transformation scene has taken place in the five years that have 
elapsed, but an encouraging fact is that one professor at any rate 
has been imported from Dublin. 

There is no lack in Cambridge of the cultivation of Quater- 
nions as an algebra, but this cultivation is not Hamiltonian, 
though an evidence of the great fecundity of Hamilton's work. 
Hamilton looked upon Quaternions as a geometrical method, and it 
is in this respect that he has as yet failed to find worthy followers 
resident in Cambridge. [The chapter contributed by Prof. Cayley 
to Prof. iTait's 3rd ed. of 'Quaternions' deals with quite a different 
subject from the rest of the treatise, a subject that deserves a 
distinctive name, say, Cayleyan Quaternions.] 

I have delayed for a considerable time the present publication 
in order at the last if possible to make it more effective. I have 
waited till I could by a more striking example than any in the 
essay shew the immense utility of Quaternions in the regions in 
which I believe them to be especially powerful. This I believe 
has been done in the *PhiL Trans.' 1892, p. 685. Certainly on 
two occasions copious extracts have been published, viz. in the 
P. R. S. R, 1890-1, p. 98, and in the ' Phil. Mag.' June 1892, p. 
477, but the reasons are simple. The first was published after the 
subject of the * Phil. Trans.' paper had been considered sufficiently 
to afford clear daylight ahead in that direction, and the second 
after that paper had actually been despatched for publication. 

At the time of writing the essay I possessed little more than 
faith in the potentiality of Quaternions, and I felt that something 
more than faith was needed to convince scientists. It was thought 
that rather than publish in driblets it were better to wait for 
a more copious shower on the principle that a well-directed heavy 
blow is more effective than a long-continued series of little 
pushes. 

Perhaps more harm has been done than by any other cause to 
the study of Quaternions in their Physical applications by a silly 
superstition with which the nurses of Cambridge are wont to 
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frighten their too timorous charges. This is the belief that the 
subject of Quaternions is difficult. It is difficult in one sense and 
in no other, in that sense in which the subject of analytical conies 
is difficult to the schoolboy, in the sense in which every subject 
is difficult whose fundamental ideas and methods are different 
from any the student has hitherto been introduced to. The only 
way to convince the nurses that Quaternions form a healthy diet 
for the young mathematician is to prove to them that they will 
" pay" in the first part of the Tripos. Of course this is an impos- 
sible task while the only questions set in the Tripos on the subject 
are in the second part and average one in two years. [This 
solitary biennial question is rarely if ever anything but an exercise 
in algebra. The very form in which candidates are invited, or at 
any rate were in my day, to study Quaternions is an insult to the 
memory of Hamilton. The monstrosity " Quaternions and other 
non-commutative algebras*' can only be paralleled by ** Cartesian 
Geometry and other commutative algebras." When I was in Cam- 
bridge it was currently reported that if an answer to a Mathe- 
matical Tripos question were couched in Hebrew the candidate 
would or would not get credit for the answer according as one 
or more of the examiners did or did not understand Hebrew, 
and that in this respect Hebrew or Quaternions were strictly 
analogous.] 

Is it hopeless to appeal to the charges ? I will try. Let me 
suppose that some budding Cambridge Mathematician has fol- 
lowed me so far. I now address myself to him. Have you ever 
felt a joy in Mathematics ? Probably you have, but it was before 
your schoolmasters had found you out and resolved to fashion you 
into an examinee. Even now you occasionally have feelings like 
the dimly remembered ones. Now and then you forget that you 
are nerving yourself for that Juggernaut the Tripos. Let me 
implore you as though your souFs salvation depended on it to let 
these trances run their utmost course in spite of solemn warnings 
from your nurse. You will in time be rewarded by a soul-thrilling 
dream whose subject is the Universe and whose organ to look 
upon the Universe withal is the sense called Quaternions. Steep 
yourself in the delirious pleasures. When you wake you will have 
forgotten the Tripos and in the fulness of time will develop into a 
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financial wreck, but in possession of the memory of that heaven- 
sent dream you will be a far happier and richer man than the 
millionest millionaire. 

To pass to earth — from the few papers I have published it 
will be evident that the subject treated of here is one I have very 
much at hearty and I think that the publication of the essay is 
likely to conduce to an acceptance of the view that it is now the 
duty of mathematical physicists to study Quaternions seriously. 
I have been told by more than one of the few who have read some 
of my papers that they prove rather stiff reading. The reasons 
for this are not in the papers I believe but in matters which have 
already been indicated. Now the present essay reproduces the 
order in which the subject was developed in my own mind. The 
less complete treatment of a subject, especially if more diffuse, is 
often easier to follow than the finished product. It is therefore 
probable that the present essay is likely to prove more easy 
reading than my other papers. 

Moreover I wish it to be studied by a class of readers who are 
not in the habit of consulting the proceedings, &c., of learned 
societies. I want the slaves of examination to be arrested and to 
read, for it is apparently to the rising generation that we must 
look to wipe off' the blot from the escutcheon of Cambridge. 

And now as to the essay itself. But one real alteration has 
been made. A passage has been suppressed in which were made 
uncomplimentary remarks concerning a certain author for what 
the writer regards as his abuse of Quaternion methods. The 
author in question would no doubt have been perfectly well able 
to take care of himself, so that perhaps there was no very good 
reason for suppressing the passage as it still represents my con- 
victions, but 1 did not want a side issue to be raised that would 
serve to distract attention from the main one. To bring the 
notation into harmony with my later papers dv and V which 
occur in the manuscript have been changed throughout to d% and 
A respectively. To facilitate printing the solidus has been freely 
introduced and the vinculum abjured. Mere slips of the pen 
have been corrected. A formal prefatory note required by the 
conditions of competition has been omitted. The Table of Con- 
tents was not prefixed to the original essay. It consists of little 
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more than a collection of the headings scattered through the 
essay. Several notes have been added, all indicated by square 
brackets and the date (1892 or 1893). Otherwise the essay 
remains absolutely unaltered. The name originally given to the 
essay is at the head of p. 1 below. The name on the title-page 
is adopted to prevent confusion of the essay with the * Phil. Mag.' 
paper referred to above. What in the peculiar caligraphy of the 
manuscript was meant for the familiar /JJ()d9 has been consistently 
rendered by the printer as ///()ds. As the mental operation of 
substituting the former for the latter is not laborious I have not 
thought it necessary to make the requisite extensive alterations in 
the proofs. 

I wish here to express my great indebtedness to Prof Tait, not 
only for having through his published works given me such 
knowledge of Quaternions as I possess but for giving me private 
encouragement at a time I sorely needed it. There was a time 
when I felt tempted to throw my convictions to the winds and 
follow the line of least resistance. To break down the solid and 
well-nigh universal scepticism as to the utility of Quaternions in 
Physics seemed too much like casting one's pearls — at least like 
crying in the wilderness. 

But though I recognise that I am fighting under Prof. Tait's 
banner, yet, as every subaltern could have conducted a campaign 
better than his general, so in some details I feel compelled to 
differ from Professor Tait. Some two or three years ago he was 
good enough to read the present essay. He somewhat severely 
criticised certain points but did not convince me on all. 

Among other things he pointed out that I sprang on the 
unsuspicious reader without due warning and explanation what 
may be considered as a peculiarity in symbolisation, I take this 
opportunity therefore of remedying the omission. In Quaternions 
on account of the non-commutative nature of multiplication we 
have not the same unlimited choice of order of the terms in a 
product as we have in ordinary algebra, and the same is true of 
certain quaternion operators. It is thus inconvenient in many 
cases to use the familiar method of indicating the connection 
between an operator and its operand by placing the former im- 
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integral of the kind. By quite elementary considerations he sees 
that while only such volume integrals as satisfy certain conditions 
are transformable into surface integrals, yet any surface integral 
which is continuous and applies to the complete boundary of any 
finite volume can be expressed as a volume integral throughout 
that volume. He is thus led to start from the surface integral 
and deduces by the briefest of processes the most general volume 
integral of the type required. Needless to say, when giving his 
demonstration he does not bare his soul in this way. He thinks 
rightly that any mathematician can at once divine the exact road 
he has followed. Where is the artificiality ? 

Let me in conclusion say that even now I scarcely dare state 
what I believe to be the proper place of Quaternions in a 
Physical education, for fear my statements be regarded as the 
uninspired babblings of a misdirected enthusiast, but I cannot 
refrain from saying that I look forward to the time when 
Quaternions will appear in every Physical text-book that assumes 
the knowledge of (say) elementary plane trigonometry. 

I am much indebted to Mr G. H. A. Wilson of Clare College, 
Cambridge, for helping me in the revision of the proofs, and take 
this opportunity of thanking him for the time and trouble he has 
devoted to the work. 

ALEX. McAULAY. 



University op Tasmania, 

HOBART. 

March 26, 1893. 
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QUATERNIONS AS A PRACTICAL INSTRU- 
MENT OF PHYSICAL RESEARCH. 



SECTION L 

Introduction. 

It is a curious phenomenon in the History of Mathematics 
that the greatest work of the greatest Mathematician of the 
century which prides itself upon being the most enlightened the 
world has yet seen, has suffered the most chilling neglect. 

The cause of this is not at first sight obvious. We have here 
little to do with the benefit provided by Quaternions to Pure 
Mathematics. The reason for the neglect here may be that 
Hamilton himself has developed the Science to such an extent as 
to make successors an impossibility. One cannot however resist 
a strong suspicion that were the subject even studied we should 
hear more from Pure Mathematicians, of Hamilton's valuable 
results. This reason at any rate cannot be assigned for the 
neglect of the Physical side of Quaternions. Hamilton has done 
but little in this field, and yet when we ask what Mathematical 
Physicists have been tempted by the bait to win easy laurels (to 
put the incentive on no higher grounds), the answer must be 
scarcely one. Prof. Tait is the grand exception to this. But well- 
known Physicist though he be, his fellow-workers for the most 
part render themselves incapable of appreciating his valuable 
services by studying the subject if at all only as dilettanti. The 
number who read a small amount in Quaternions is by no means 
small, but those who get further than what is recommended by 
Maxwell as imperatively necessary are but a small percentage of 
the whole. 

I cannot help thinking that this state of affairs is owing chiefly 
to a prejudice. This prejudice is well seen in Maxwell's well- 
known statement — " I am convinced that the introduction of the 

M. 1 



2 INTRODUCTION. 

ideas, as distinguished from the operations and methods of Quater- 
nions, will be of great use to us in all parts of our subject/'* 
Now what I hold and what the main object of this essay is to 
prove is that the "operations and methods" of Quaternions are as 
much better qualified to deal with Physics than the ordinary ones 
as are the *' ideas". 

But, what has produced this notion, that the subject of 
Quaternions is only a pretty toy that has nothing to do with the 
serious work of practical Physics? It must be the fact that it has 
hitherto produced few results that appeal strongly to Physicists. 
This I acknowledge, but that the deduction is correct I strongly 
disbelieve. As well might an instrument of which nobody has 
attempted to master the principles be blamed for not being of 
much use. Workers naturally find themselves while still inex- 
perienced in the use of Quaternions incapable of clearly thinking 
through them and of making them do the work of Cartesian Geo- 
metry, and they conclude that Quaternions do not provide suitable 
treatment for what they have in hand. The fact is that the sub- 
ject requires a slight development in order readily to apply to the 
practical 'Consideration of most physical subjects. The first steps 
of this, which consist chiefly in the invention of new symbols of 
operation and a slight examination of their chief properties, I 
have endeavoured to give in the following pages. 

I may now state what I hold to be the mission of Quaternions 
to Physics. I believe that Physics would advance with both 
more rapid and surer strides were Quaternions introduced to 
serious study to the almost total exclusion of Cartesian Geo- 
metry, except in an insignificant way as a particular case of 
the former. All the geometrical processes occurring in Physical 
theories and general Physical problems are much more grace- 
ful in their Quaternion than in their Cartesian garb. To illus- 
trate what is here meant by "theory" and "general problem" 
let us take the case of Elasticity treated below. That by the 
methods advocated not only are the already well-known results of 
the general theory of Elasticity better proved, but more general 
results are obtained, will I think be acknowledged after a perusal 
of §§ 12 to 21 below. That Quaternions are superior to Cartesian 
Geometry in considering the general problems of (1) an infinite 
isotropic solid, (2) the torsion and bending of prisms and cylinders, 
(3) the general theory of wires, I have endeavoured to shew in 

* ElecU and Map. Vol. i. § 10. 
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§§ 22 — 33. But for particular' problems such as the torsion 
problem for a cylinder of given shape, we require of course the 
various theories specially constructed for the solution of particular 
problems such as Fourier's theories, complex variables, spherical 
harmonics, &c. It will thus be seen that I do not propose to 
banish these theories but merely Cartesian Geometry. 

So mistaken are the common notions concerning the preten- 
sions of advocates of Quaternions that I was asked by one well- 
known Mathematician whether Quaternions furnished methods for 
the solution of differential equations, as he asserted that this was 
all that remained for Mathematics in the domain of Physics ! 
Quaternions can no more solve differential equations than Carte- 
sian Geometry, but the solution of such equations can be performed 
as readily, in fact generally more so, in the Quaternion shape as in 
the Cartesian. But that the sole work of Physical Mathematics 
to-day is the solution of differential equations I beg to question. 
There are many and important Physical questions and extensions 
of Physical theories that have little or nothing to do with such 
solutions. As witness I may call attention to the new Physical 
work which occurs below. 

If only on account of the extreme simplicity of Quaternion 
notation, large advances in the parts of Physics now indicated, are 
to be expected. Expressions which are far too cumbrous to be of 
much use in the Cartesian shape become so simple when translated 
into Quaternions, that they admit of easy interpretation and, what 
is perhaps of more importance, of easy manipulation. Compare for 
instance the particular case of equation (15 m) § 16 below when 
P = with the same thing as considered in Thomson and Tait's 
Nat Phil., App. C. The Quaternion equation is 

p\8Vi^(jwA = 0. 

The Cartesian exact equivalent consists of Thomson and Tait's 
equations (7), viz. 

dw doL dw da] 
db dz dc dy 

d {^dw doL dw da dw (da 




dy 

d^ 
dz 



dB dy da dz dc \dx 

dw da dw da dw (da ^ > r _ ^ 
dC dz da dy db \dx ' 



and two similar equations. 
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4 INTRODUCTION. 

Many of the equations indeed in the part of the essay where 
this occurs, although quite simple enough to be thoroughly useful 
in their present form, lead to much more complicated equations 
than those just given when translated into Cartesian notation. 

It will thus be seen that there are two statements to make 
good : — (1) that Quaternions are in such a stage of development 
as already to justify the practically complete banishment of Carte- 
sian Geometry from Physical questions of a general nature, and 
(2) that Quaternions will in Physics produce many iiei(; results that 
cannot be produced by the rival and older theory. 

To establish completely the first of these propositions it would 
be necessary to go over all the ground covered by Mathematical 
Physical Theories, by means of our present subject, and compare 
the proofs with the ordinary ones. This of course is impossible in 
an essay. It would require a treatise of no small dimensions. 
But the principle can be followed to a small extent. I have there- 
fore taken three typical theories and applied Quaternions to most 
of the general propositions in them. The three subjects are those 
of Elastic Solids, with the thermodynamic considerations necessary, 
Electricity and Magnetism, and Hydrodynamics. It is impossible 
witliout greatly exceeding due limits of space to consider in addi- 
tion, Conduction of Heat, Acoustics, Physical Optics, and the 
Kinetic Theory of Gases. With the exception of the first of these 
subjects I do not profess even to have attempted hitherto the 
desired applications, but one would seem almost justified in argu- 
ing that, since Quaternions have been found so applicable to the 
subjects considered, they are very likely to prove useful to about 
the same extent in similar theories. Again, only in one of the 
subjects chosen, viz.. Hydrodynamics, have I given the whole of 
the general theory which usually appears in text-books. For 
instance, in Electricity and Magnetism I have not considered Elec- 
tric Conduction in three dimensions which, as Maxwell remarks, 
lends itself very readily to Quaternion treatment, nor Magnetic 
Induction, nor the Electro-Magnetic Theory of Light. Again, I 
have left out the consideration of Poynting's theories of Elec- 
tricity which are very beautifully treated by Quaternions, and I 
felt much tempted to introduce some considerations in connection 
with the Molecular Current theory of Magnetism. With similar 
reluctance I have been compelled to omit many applications in 
the Theory of Elastic Solids, but the already too large size of the 
essay admitted of uo other course. Notwithstanding these omis- 
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sions, I think that what I liave done in this part will go far to 
bear out the truth of the first proposition I have stated above. 

But it is the second that I would especially lay stress upon. 
In the first it is merely stated that Cartesian Geometry is an anti- 
quated machine that ought to be thrown aside to make room for 
modern improvements. But the second asserts that the improved 
machinery will not only do the work of the old better, but will 
also do much work that the old is quite incapable of doing at all. 
Should this be satisfactorily established and should Physicists in 
that case still refuse to have anything to do with Quaternions, 
they would place themselves in the position of the traditional 
workmen who so strongly objected to the introduction of ma- 
chinery to supplant manual labour. 

But in a few months and synchronously with the work I have 
already described, to arrive at a large number of new results is too 
much to expect even from such a subject as that now under dis- 
cussion. There are however some few such results to shew. I 
have endeavoured to advance each of the theories chosen in at 
least one direction. In the subject of Elastic Solids I have ex- 
pressed the stress in terms of the strain in the most general case, 
i.e. where the strain is not small, where the ordinary assumption 
of no stress-couple is not made and where no assumption is made 
as to homogeneity, isotropy, &c. I have also obtained the equations 
of motion when there is given an external force and couple per 
unit volume of the unstrained solid. These two problems, as will 
be seen, are by no means identical. In Electrostatics I have con- 
sidered the most general mechanical results flowing from Maxwell's 
theory, and their explanation by stress in the dielectric. These 
results are not known, as might be inferred from this mode of 
statement, for to solve the problem we require to know forty-two 
independent constants to express the properties of the dielectric at 
a given state of strain at each point. These are the six coefficients 
of specific inductive capacity and their thirty-six differential coeffi- 
cients with regard to the six coordinates of pure strain. But, as 
far as I am aware, only such particular cases of this have already 
been considered as make the forty-two constants reduce at most to 
three. In Hydrodynamics I have endeavoured to deduce certain 
general phenomena which would be exhibited by vortex-atoms 
acting upon one another. This has been done by examination of 
an equation which has not, I believe, been hitherto given. The 
result of this part of the essay is to lead to a presumption against 
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Sir William Thomson's Vortex- Atom Theory and in favour of 
Hicks's. 

As one of the objects of this introduction is to give a bird's-eye 
view of the merits of Quaternions as opposed to Cartesian Geome- 
try, it will not be out of place to give side by side the Quaternion 
and the Cartesian fonns of most of the new results I have been 
speaking about. It must be premised, as already hinted, that the 
usefulness of these results must be judged not by the Cartesian 
but by the Quaternion form. 



Elasticity. 

Let the point {x, y, z) of an elastic solid be displaced to 
{x\ y\ z). The strain at any point that is caused may be supposed 
due to a pure strain followed by a rotation. In Section III. below, 
this pure strain is called ^^^ Let its coordinates be e, f, g, 
a/2, 6/2, c/2; i.e. if the vector (^, 77, ^) becomes (^', 77', f) by means 
of the pure strain, then 

f = e? + ic77 + 46?, 
&c., &c. 

Thus when the strain is small Cyfyg reduce to Thomson and Tait's 
1 + e, 1 +/, l+g and a, b, c are the same both in their case and 
the present one. Now let the coordinates of ^, § 16 below, be 
E, F, 0, A 12, i^/2, 0/2. Equation (15), § 16 below, viz. 

gives in our present notation 

E = e' + cV4 + 6V4 = (cLjc'Id^y + (dy'/dxf + (d//cfo)^ 

&c., &c. 

A=.a(f+g)-{- 6c/2 = 2 [(dx'/dy) {dxjdz) + {dy'jdy) (dy/dz) 

+ (dz'/dy){dz'/dz)], 
&c., &c. 

which shew that the present E, F, (?, A/ 2, 5/2, (7/2 are the 
Ay Bj G, a, 6, c of Thomson and Tait's i\^a^. Phil., App. C. 

* This result is one of Tait's (Quaternions § 365 where he has <p'(p=m^). It is 
given here for completeness. 
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J ill) ='^"' ^^ ^''•' 

J (^'^') = J.r J (f^) = J.V &^. &C., &C., &C. 

I have shewn in § 14 below that the stress-couple is quite 
independent of the strain. Thus we may consider the stress to 
consist of two parts — an ordinary stress PQRSTU as in Thomson 
and Tait's Nat, Phil, and a stress which causes a couple per unit 
volume L'M'N\ The former only of these will depend on strain. 
The result of the two will be to cause a force (as indeed can be 
seen from the expressions in § 13 below) per unit area on the x- 
interface P, TJ^-N'\% T-Mj^, and so for the other interfaces. If 
Zr, M, N be the external couple per unit volume of the unstrained 
solid we shall have ' 

for the external couple and the stress- couple are always equal and 
opposite. Thus the force on the a;-interface becomes 

P, U-NI2J, T + M/2J 

and similarly for the other interfaces. 

To express the part of the stress (P &c.) which depends on 
the strain in terms of that strain, consider w the potential energy 
per unit volume of the unstrained solid as a function of E &c. 
In the general thermodynamic case w may be defined by saying 
that w X (the element of volume) = (the intrinsic energy of the 
element) — (the entropy of the element x its absolute temperature 
X Joule's coefficient). Of course w may be, and indeed is in 
§§ 14, 15 below, regarded as a function of e &c. 

The equation for stress is (15 6) § 1(5 below, viz., 

The second of the expressions is in terms of the strain and the 
third in terms of the displacement and its derivatives. In our 
present notation this last is 
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dy 

^ dx' djx' dw ^ dx' dx' dw ^ dx' dx dw 
dy dz dA dz dx dB dx dy dC 

&c., &c. 



JS _ dy^ dz' dw d'lf dz' dw dy' dz* dw 
~2~'^dx dx dE'^ dy dy dF^ dz dz dG 

(dy' dz dy' dz\ dw /dy' dz' dy' dz 
\dy dz dz dyJdA \dz dx dx dz 



(dy' dz' dy' dz'\ dw 
\dx dy dy dx) dC 



dy dy 

&c., &c. 

In § 14 I also obtain this part of the stress explicitly in terms 
of Cy f, g, a, 6, c, of i(; as a function of these quantities and of the 
axis and amount of rotation. But these results are so very com- 
plicated in their Cartesian shape that it is quite useless to give 
them. 

To put down the equations of motion let X^, Y^, Z^ be the 
force due to stress on what before strain was unit area per- 
pendicular to the axis of x. Similarly for Xy, &c. Next suppose 
that X, F, Z is the external force per unit volume of the un- 
strained solid and let D be the original density of the solid. Then 
the equation of motion (15 ?i) § 16 a below, viz. 

D/3' = P + tA, 

gives in our present notation 

X + dXJdx + dXJdy + dXJdz = xDy 

&c., &c. 

It remains to express X^ &c. in terms of the displacement and 
LMN, This is done in equation (15 i) § 16 below, viz. 

In our present notation this consists of the following nine 
equations : 

* The second term on the right contains in full the nine terms corresponding to 
(Jy2N-Ji^M)j2J. Quaternion notation is therefore here, as in nearly all cases 
which occur in Physics, considerably more compact even than the notations of 
determinants or Jacobians. 
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'"''' \dEdx '^ dCdy '^dBdzJ'^ 2J 

y __a /^^ ^ , ^^ ^^-s:' dw d/\ JuM — J12L 
'~ [dEd^ ^dCd^ '^dBdzJ'^ 2J ' 

and six similar equations. 

We thus see that in the case where LMN are zero, our present 
X'^t Xyy X, are the PQR of Thomson and Tait's Nat Phil, App. C 
(d), and therefore equations (7) of that article agree with our 
equations of motion when we put both the external force and the 
acceleration zero. 

These are some of the new results in Elasticity, but, as I have 
hinted, there are others in §§ 14, 15 which it would be waste of 
time to give in their Cartesian form. 



Electricity, 

In Section IV. below I have considered, as already stated, the 
most general mechanical results flowing from Maxwell's theory of 
Electrostatics. I have shewn that here, as in the particular cases 
considered by others, the forces, whether per unit volume or per 
unit surface, can be explained by a stress in the dielectric. It is 
easiest to describe these forces by means of the stress. 

Let the coordinates of the stress be PQRSTIT, Then F^F^F^ 
the mechanical force, due to the field per unit volume, exerted 
upon the dielectric where there is no discontinuity in the stress, is 
given by 

F^ = dP/dx + dU/dy + dTjdz, &c., &c. 

and (i, m, n) being the direction cosines of the normal to any 
surface, pointing away from the region considered 

F; =^^[lP + mU+ nT], - [ ]„ &c., &c., 

where a, b indicate the two sides of the surface and -F/, i^Y, F^ is 
the force due to the field per unit surface. 

It remains to find P &c. Let X, Y, Z be the electro-motive 
force, a, yS, 7 the displacement, iv the potential energy per unit 
volume and K^^, K^y, K^, Ky^, K^y Kj,y the coefficients of specific 
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inductive capacity. Let 1 -H e, 1 +/, 1+5^, a/2, 6/2, c/2 denote the 
pure part of the strain of the medium. The K's will then be 
functions of e &c. and we must suppose these functions known, or 
at any rate we must assume the knowledge of both the values of 
the K*8 and their differential coefficients at the particular state of 
strain in which the medium is when under consideration. The 
relations between the above quantities are 

i7roL = K^ + K^Y + K,^, &c., &c. 

w = (Xa + Y^ + Zy)l2 
= (K^^X' + K^Y' + Kj;' + 2K,JZ + 2KJX 

+ 2K^^XY)/S7r, 

It is the second of these expressions for w which is assumed 
below, and the differentiations of course refer only to the K's, 
The equation expressing P &c. in terms of the field is (21) § 40 
below, viz. 

cjxo = — ^ FDeaE — ^(JwcOy 

which in our present notation gives the following six equations 

P= - 1( _ aX + ^F+ 7^ - dw/de, &c., &c., 
S==i(^Z+yY)- dwjda, &c., &c. 

I have shewn in §§ 41 — 45 below that these results agree with 
particular results obtained by others. 

Hydrodynamics, 

The new work in this subject is given in Section VI. — " The 
Vortex- Atom Theory." It is quite unnecessary to translate the 
various expressions there used into the Cartesian form. I give 
here only the principal equation in its two chief forms, equation 
(9) § 89 and equation (11) § 90, viz. 

P+v- 0-72 + (47r)-YJ/(/Sf(TTV« + lidmldt) d% = H, 

F + v- a' 12 + (47r)-VJ/ {dsS^u ( Var - ma) + ud {md%)ldt} = H, 

In Cartesian notation these are 
Jdplp+V + q'l2 

- (47r)"V//{2 [(x' - x) (tvv - y?) + . . . +. . .]/^'' + {dcldt)lr]dxdy'dz = H. 

jdplp+V^-ifl2 

- (47r)"7JJ [{x - x) [2 (i^T; - ^;f ) - c^O + . . . + . . . }/»'' • dx'dydz 

- (47r)-V// [d (cdxdyW)ldt]lr = K 
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The fluid here considered is one whose motion is continuous 
from point to point and which extends to infinity. The volume 
integral extends throughout space. The notation is as usual. 
It is only necessary to say that H is a, function of the time only, 
r is the distance between the points a?', y\ z! and x,y, z\ 

c = dujdx + dvjdy 4- dwjdz ; 

djdt is put for diflferentiation which follows a particle of the fluid, 
and djdt for that which refers to a fixed point. 

The explanation of the unusual length of this essay, which I 
feel is called for, is contained in the foregoing description of its 
objects. If the objects be justifiable, so must also be the length 
which is a necessary outcome of those objects. 



SECTION II. 

Quaternion Theorems. 
Definitions. 

1. As there are two or three symbols and terms which will be 
in constant use in the following pages that are new or more general 
in their signification than is usual, it is necessary to be perhaps 
somewhat tediously minute in a few preliminary definitions and 
explanations. 

A function of a variable in the following essay is to be under- 
stood to mean anything which depends on the variable and which 
can be subjected to mathematical operations; the variable itself 
being anything capable of being represented by a mathematical 
symbol. In Cartesian Geometry the variable is generally a single 
scalar. In Quaternions on the other hand a general quaternion 
variable is not infrequent, a variable which requires 4 scalars for 
its specification, and similarly for the function. In both, however, 
either the variable or the function may be a mere symbol of 
operation. In the following essay we shall frequently have to 
speak of variables and functions which are neither quaternions nor 
mere symbols of operation. For instance K in § 40 below requires 
6 scalars to specify it, and it is a function of -^/r which requires 
6 scalars and p which requires 3 scalars. When in future the 
expression " any function " is used it is always to be understood in 
the general sense just explained. 

We shall frequently have to deal with functions of many inde- 
pendent vectors, and especially with functions which are linear in 
each of the constituent vectors. These functions merely require 
to be noticed but not defined. 

Hamilton has defined the meaning of the symbolic vector V 
thus : — 

dx ^ dy dz' 
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where i, j, k are unit vectors in the directions of the mutually 
perpendicular axes x, y, z, I have found it necessary somewhat to 
expand the meaning of this symbol. When a numerical suffix 
1, 2,... is attached to a V in any expression it is to indicate that the 
differentiations implied in the V are to refer to and only to other 
symbols in the same expression which have the same suffix. After 
the implied differentiations have been performed the suffixes are 
of course removed. Thus Q (a, )8, 7, S) being a quaternion func- 
tion of any four vectors a, )8, 7, S, linear in each 

and again 

It is convenient to reserve the symbol A for a special meaning. 
It is to be regarded as a particular form of V, but its differentia- 
tions are to refer to all the variables in the term in which it 
appears. Thus Q being as before 

If in a linear expression or function Vj and p^ (p being as usual 
= IX -\-jy + kz) occur once each they can be interchanged. Simi- 
larly for Vg and p^. So often does this occur that I have thought 
it advisable to use a separate symbol f, for each of the two Vj and 
Pj, fg for each of the two V^ and p^ and so for fg, &c. If only one 
such pair occur there is of course no need for the suffix attached to 
f. Thus f may be looked upon as a symbolic vector or as a single 
term put down instead of three. For Q (a, /8) being linear in each 
of the vectors a, /8 

Q (r. r) = Q (V., p.) = Q a i) + q uj) + q {k k) (i). 

There is one more extension of the meaning of V to be given. 

* These meanings for Vi» V2--^ ^ "^^^ i^ ^ paper on '* Some General Theorems 
in Quaternion Integration," in the Mess, of Math. Vol. xiv. (18»4), p. 26. The 
investigations there given are for the most part incorporated below. [Note added, 
1892, see preface as to the alteration of y' into A.] 
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u, V, w being the rectangular coordinates of any vector o*, ^V is 
defined by the equation 

— _ . d . d , d 

"^ du *' dv dw' 

To <,V of course are to be attached, when necessary, the suffixes 
above explained in connection with V. Moreover just as for 
Vj, p^ we may put f, f so also for <^Vj, a^ may we put the same. 

With these meanings one important result follows at once. 
The Vj s, Vj's, cfcc, obey all the laws of ordinary vectors whether 
vnth regard to multiplication or addition, for the coordinates 
didx, djdy, d/dz of any V obey with the coordinates of any vector 
or any other V all the laws of common algebra. 

Just as <^V may be defined as a symbolic vector whose coordi- 
nates are d/du, d/dv, d/dw so <f> being a linear vector function of 
any vector whose coordinates are 

{afi^c^ ajb^c^ ajb^c^) (i.e. ^i = a,i + bj + cjc, &c.). 
^Q* is defined as a symbolic linear vector function whose coordi- 
nates are 

(d/da^y d/d\y d/dc^, djda^, djdb^y d/dc^y d/da^y d/d\y d/dc^), 

and to ^(j is to be applied exactly the same system of suffixes as 
in the case of V. Thus q being any quaternion function of ^, and 
6) any vector 

4>Gi^ •?! = "" {i(iqlda^ +jdq/db^ + kdq/dc^) Siay 

- {idqjda^ +jdq/db^ + kdqjdc^) Sj(o 

— (idq/da^ -\-jdqldb^ -h kdq/dc^) Skoo. 

The same symbol ^(j is used without any inconvenience with a 
slightly different meaning. If the independent variable ^ be a 
self-conjugate linear vector function it has only six coordinates. 
If these are PQRSTU (i.e. 4>i = Pi + Uj + Tk, &c.) ^a is defined 
as a self-conjugate linear vector function whose coordinates are 

{djdP, d/dQy d/dRy i^djdSy \d j dT y ^d/dU). 

We shall frequently have to compare volume integrals with 
integrals taken over the bounding surface of the volume, and 
again surface integrals with integrals taken round the boundary 
of the surface. For this purpose we shall use the following nota- 
tions for linear, surface and volume integrals respectively JQdp, 

* I have used an inverted D to indicate the analogy to Hamilton's inverted A. 
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JjQd^y fJfQd% where Q is any function of the position of a point. 
Here dp is a vector element of the curve, dS a vector element of 
the surface, and d% an element of volume. When comparisons 
between line and surface integrals are made we take dS in such a 
direction that dp is in the direction of positive rotation round the 
element dX close to it. When comparisons between surface and 
volume integrals are made dS is always taken in the direction 
away from the volume which it bounds. 



Properties of f. 

2. The property of f on which nearly all its usefulness 
depends is that if <r be any vector 

cr = - ?/S?(r, 

which is given at once by equation (1) of last section. 

This gives a useful expression for the conjugate of a linear 
vector function of a vector. Let <^ be the function and 6), t any 
two vectors. Then ^' denoting as usual the conjugate of ^ we 
have 

S(0<f>T = ST(f>'(0, 

whence putting on the left t = — fSffr we have 

St (- ^S(o(l>0 = St^'6), 

or since t is quite arbitrary 

fft) = -fSft>(^? (2). 

From this we at once deduce expressions for the pure part ^co 
and the rotational part Veo) of <f>(o by putting 

{<l> -<!>') CO = -(i>^Scoi;-\-^Sco(f>^=vv^(f>^. CO = 2V€coy"^^' 

And all the other well-known relations between cf) and <f> are at 
once given e.g. S^<f>^ = S^^' ^, i.e. the *' convergence " of ^ = the 
" convergence " of <f>\ 

3. Let Q (\, fji) be any function of two vectors which is linear 
in each. Then if cfyco be any linear vector function of a vector co 
given by 

we have Q(f, (^0 = SQ(- C^^ra,/3) = S(3(a, )8) ) ^ ^' 
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or more generally 

<2 (?. ^X?) = 2(2 (x'a, )8) (4a). 

To prove, it is only necessary to observe that 

<^X?= - 2/3Sax?= - 2/9S?x'«' 
and that — ^S^x'a — x'«- 

As a particular case of eq. (4) let ^ have the self-conjugate value 

Then Q(?,<^?) = i2{Q(a,/3) + Q(Aa)}j ^ ^' 

or if Q (\, /j) is symmetrical in \ and fi 

(2(r,<^?) = SQ(a,/3) (6). 

The application we shall frequently make of this is to the case 
when for a we put Vj and for y8, cr,, where a is any vector function 
of the position of a point. In this case the first expression for (f) 
is the strain function and the second expression the pure strain 
function resulting from a small displacement a at every point. 
As a simple particular case put Q (\, fi) = S\/jl so that Q is sym- 
metrical in \ and /x. Thus (f> being either of these functions 

Another important equation is 

*Q(?.«^r) = <3(f?,0 (6a). 

This is quite independent of the form of (f>. To prove, observe 
that by equation (2) 

and that — fSft^f, = <f>^^. Thus we get rid of f and may now 
drop the suffix of fj and so get eq. (Qa), [Notice that by means 
of (6a), (4a) may be deduced from (4); for by (6a) 

Q (?. <t>xO = Q (x'f. </>?) = SQ (x'«, /8) by (4).] 

3a. A more important result is the expression for <fr^cD in 
terms of <t>. We assume that 

S(f>\(l>fi<l>v = mSXfiVy 
where X, fi, v are any three vectors and m is a scalar independent 

* [Note added, 1892. For practice it is convenient to remember this in words : — 
A term in which f and <f>^ occur is unaltered in value by chanping them into 0'f 
and ^respectively.] 
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of these vectors. Substituting f^, f^, fg for \, /a, z; and multiplying 

'Sf^.^a'S*^?.^?;*/'?, = 6«i (66), 

which gives m in terms of (f>. That fi^?,§'2&'^?i?2?8 = 6 is seen by 
getting rid of each pair of f*s in succession thus : — 

= - 2?,?, = 6. 
Next observe that 

Multiplying by Ff^f^ *^^ again on the right getting rid of the 
^s we have 

Fr.r^<^a><^(;K, = -2m« (6c), 

whence from equation (66) 
or changing « into <^-'a, 

By equation (6a) of last section we can also put this in the form 

^ ''—muMKMM. ^ ^' 

so that ifr^fo is obtained explicitly in terras of <^ or <f>. 

Equation {'^c) or (6d) can be put in another useful form which 
is more analogous to the ordinary cubic and can be easily deduced 
therefrom, or *less easily from (Gd), viz. 

= (6/). 

* [Note added, 1892. The cubic may be obtained in a more useful form from 
the equation wStM^S<t>tM^<t>^^= -^V^i^^S(P(^<pt,(f>^^ thus 

Again <pti^ . ipu)(t> ^^^l: 1^2= Hi^ - <P<^V . <i>^^Vt^t2=<t>^\^<t>^ ( " ti^t-M^-^ ^^^^M^) 

= - (ri'S'ri0cu) s^^4>^.+ct> (fi-sri^ro) ^0^ 

Hence ^^w - m"<l>^(a + m'^w - wiw = 0, 

where 6m = S^^ ^^iS4>^^4>^^4>^^ 

2m' = - SVi:,t2V4>^M2 
m"=-Smi 
M. 2 
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As a useful particular case of equation (6d) we may notice 
that by equation (4) § 3 if 

<f)Q) = — o-jiSwVj 
<^-'a) = - 3 FV,V^a,o-,o-,/<SrV.V.V^<7,<7.<7, (Qg), 

and if>'-% = -3F<r,o-,/Sft)V^V^SV,V^V,So-,<r,<r, (Gh). 

4. Let ^, ■^ be two linear vector functions of a vector. 
Then if 

where % is a quite arbitrary linear vector function 

for we may put ;j^f =Tfiifft) where t and co are arbitrary vectors, so 

that ' 

ST(f)Q) = SryfrG), 

or ^(o = ylr(o. 

Similarly* if ^ and yjr are both self-conjugate and % is a quite 
arbitrary self-conjugate linear vector function the same relation 
holds as can be seen by putting 



Fundamental property of q. 

5. Just as the fundamental property of <^V is that, Q being 
any function of a 

so we have a similar property of q. Q being any function of <^ a 
linear vector function 

SQ = - Q.'SfS<^?«a.? (7). 

The property is proved in the same way as for V, viz. by ex- 
panding SS(f>^(j^^ in terms of the coordinates of <^q. First let <f> 
be not self-conjugate, and let its nine coordinates be 

{aj},c^aj)^c^aj)^c^). 

* [Note added, 1892. The following slightly more general statement is a practi- 
cally much more convenient form of enunciation : if Sxt<P^= >^xi^^^f where x is a 
perfectly arbitrary self -con jugate and <p and yf/ are not necessarily self-conjugate 
then ^=1?]. 
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Thus 

= Sa^dQ/da, + Bb^dQ/db^ + Sc^dQ/dc^, 
+ Sa,dQlda, + Sb,dQ/db, + ^,dQ/dc,, 
+ Ba^dQ/da^ + Sb,dQ/db^ + Sc^dQ/dc, = 8Q. 

The proposition is exactly similarly proved when ^ is self-con- 
jugate. 



Theorems in Integration, 

6. Referring back to § 1 above for our notation for linear 
surface and volume integrals we will now prove that if Q be any 
linear function of a vector * 

fQdp = fJQ (Vd%^) (8), 

fJQdX^mMs (9). 

To prove the first divide the surface up into a series of 
elementary parallelograms by two families of lines — one or more 
members of one family coinciding with the given boundary, — 
apply the line integral to the boundary of each parallelogram and 
sum for the whole. The result will be the linear integral given in 
equation (8). Let the sides of one such parallelogram taken in 
order in the positive direction be or, /8 -I- /8', — a — a', — y8 ; so that 
a and y8' are infinitely small compared with a and )8, and we have 
the identical relation 

The terms contributed to JQdp by the sides a and — a — a! will 
be (neglecting terms of the third and higher orders of small 
quantities) 

Qoi-Qa- Qol + Q,ctfif^V, = - Qa + Q,aS^^ ,. 

Similarly the terms given by the other two sides will be 

* These two propositionR are generalisations of what Tait and Hicks have from 
time to time proved. They were first given in the present form by me in the 
article already referred to in § 1 above. In that paper the necessary references are 
given. 

2—2 
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SO that remembering that ^8' — a' = and therefore Q^ — Qa! = 
we have for the whole boundary of the parallelogram 

where d% is put for FayS. Adding for the whole surface we get 
equation (8). 

Equation (9) is proved in an exactly similar way by splitting 
the volume up into elementary parallelepipeda by three families 
of surfaces one or more members of one of the families coinciding 
with the given boundary. If or, ^, 7 be the vector edges of one 
such parallelepiped we get a term corresponding to Qff — Qol viz. 

Q (vector sum of surface of parallelepiped) = 0, 

and we get the sum of three terms corresponding to 

above, viz. 

- Q. ( F/37) fifaV. - Q, ( Vri) -Sf/SV, - Q, ( Va^) SyV^ = - Qy,8a^y, 

whence putting Sa/Sy = — ds we get equation (9). 

7. It will be observed that the above theorems have been 
proved only for cases where we can put dQ = — Qi SdpV^ i.e. when 
the space fluxes of Q are finite. If at any isolated point they are 
not finite this point must be shut off from the rest of the space by 
a small closed surface or curve as the case may be and this surface 
or curve must be reckoned as part of the boundary of the space. 
If at a surface (or curve) Q has a discontinuous value so that its 
derivatives are there infinite whereas on each side they are 
finite, this surface (or curve) must be considered as part of the 
boundary and each element of it will occur twice, i.e. once for the 
part of the space on each side. 

In the case of the isolated points, if the surface integral or line 
integral round this added boundary vanish, we can of course cease 
to consider these points as singular. Suppose Q becomes infinite 
at the point p = a. Draw a small sphere of radius a and also a 
sphere of unit radius with the point a for centre, and consider the 
small sphere to be the added boundary. Let dS' be the element 
of the unit sphere cut off by the cone which has a for vertex and 
the element dS of the small sphere for base. Then d% = a'dS' 
and we get for the part of the surface integral considered a^fJQdid^' 
where Qax is the value of Q at the element dS. If then 
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the point may be regarded as not singular. If the limiting ex- 
pression is finite the added surface integral will be finite. If the 
expression is infinite the added surface integral will be generally 
but not always infinite. Similarly in the case of an added line 
integral if Li t {p-a)=o T (p — a) QU (p — a) is zero or finite, the added 
line integral will be zero or finite respectively (of course including 
in the term finite a possibility of zero value). If this expression 
be infinite, the added line integral will generally also be infinite. 

This leads to the consideration of potentials which is given 
in §9. 

8. Some particular cases of equations (8) and (9) which 
(except the last) have been proved by Tait, are very useful. 
First put Q = a simple scalar P. Thus 

JPdp=JJVdXVP (10), 

JfPd%=JJJVPd% (11). 

If P be the pressure in a fluid— //PdS is the force resulting from 
the pressure on any portion and equation (11) shews that — VP is 
the force per unit volume due to the same cause. Next put 
Qco = Scoa and Vcoa. Thus 

JSdpa=JSSdl,V<T (12), 

fVdpa- =JJViVd%V.<T)=JJdXSWa-iiV^Sdl.<r^... (IS), 

fJSd^a=ffJ8Vad% (14), 

fJVdt<r = JjJWadS ....(15). 

Equations (12) and (14) are well-known theorems, and (13) and 
(15) will receive applications in the following pages. Green's 
Theorem with Thomson's extension of it are, as indeed has been 
pointed out by Tait particular cases of these equations. 

Equations (14) and (15) applied to an element give the well- 
known physical meanings for S^a and V^a*. The first is 
obtained by applying (14) to any element, and the second (regard- 
ing (7 as a velocity) is obtained by applying (15) to the element 
contained by the following six planes each passing infinitely near 
to the point considered — (1) two planes containing the instanta- 
neous axis of rotation, (2) two planes at right angles to this axis, 
and (3) two planes at right angles to these four. 



* 



[Note added, 1892. Let me disarm criticism by confessing that what follows 
concerning Vya is nonsense.] 
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One very frequent application of equation (9) may be put in 
the following form: — Q being any linear function (varying from 
point to point) of R^ and V^, R being a function of the position of 
a point 

m (Ri, ^.) d% = - ///Q. {R, V,) d% + m (B, dt) . . .(16). 



Potentials. 

9. We proceed at once to the application of these theorems 
in integration to Potentials. Although the results about to be 
obtained are well-known ones in Cartesian Geometry or are easily 
deduced from such results it is well to give this quaternion method 
if only for the collateral considerations which on account of their 
many applications in what follows it is expedient to place in this 
preliminary section. 

If R is some function oi p — p where p is the vector coordinate 
of some point under consideration and p the vector coordinate of 
any point in space, we have 

Now let Q{R) be any function of R, the coordinates of Q being 
functions of p only. Consider the integral ///Q (jB) d% the variable 
of integration being p (p being a constant so far as the integral is 
concerned). It does not matter whether the integral is a volume, 
surface or linear one but for conciseness let us take it as a volume 
integral. Thus we have 

P'V///Q (B) dS = ///p'VQ (B) ds = - ///, V,Q (B,) d%. 

Now pV operating on the whole integral has no meaning so we 
may drop the affix to the V outside and always understand p\ 
Under the integral sign however we must retain the affix p or p' 
unless a convention be adopted. It is convenient to adopt such a 
convention and since Q will probably contain some pV but cannot 
possibly contain a pV we must assume that when V appears with- 
out an affix under the integral sign the affix p is understood. 
With this understanding we see that when V crosses the integral 
sign it must be made to change sign and refer only to the part we 
have called R. Thus 

V///Q (B) d8 = - fJJ^.Q (E.) d% (17). 
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Generally speaking R can and will be put as a function of 
T~^ {p — p) and for this we adopt the single symbol u. Both this 
symbol and the convention just explained will be constantly re- 
quired in all the applications which follow. 

10. Now let Q be any function of the position of a point. 
Then the potential of the volume distribution of Q, say q, is given 

q = jjjuQd% (18), 

the extent of the volume included being supposed given. We 
may now prove the following two important propositions 

V'q = 4^wQ (19), 

JJSdXV.q=^7rJJJQd^ (20). 

The latter is a corollary of the former as is seen from equation (9) 
§ 6 above. 

Equation (19) may be proved thus. If P be any function of 
the position of a point which is finite but not necessarily continu- 
ous for all points 

is always fiuite and if the volume over which the integral extends 
is indefinitely diminished, so also is the expression now under 
consideration, and this for the point at which is this remnant of 
volume. This in itself is an important proposition. The expres- 
sion, by equation (17), = — JJJ^uPd% and both statements are 
obviously true except for the point p. For this point we have 
merely to shew that the part of the volume integral just given 
contributed by the volume indefinitely near to p vanishes with 
this volume. Divide this near volume up into a series of element- 
ary cones with p' for vertex. If r is the (small) height and day the 
solid vertical angle of one of these cones, the part contributed by 
this cone is approximately U(p — p)Pprdco/S where Pp' is the 
value of P at the point p. The proposition is now obvious. 

Now since 

we see that the only part of the volume integral JffuQd% which 
need be considered is that given by the volume in the immediate 
neighbourhood of p, for at all points except p, V^u = 0. Consider 
then our volume and surface integrals only to refer to a small 
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sphere with p' for centre and so small that no point is included at 
which Q is discontinuous and therefore VQ infinite. (This last 
assumes that Q is not discontinuous at p. In the case when Q is 
discontinuous at p no definite meaning can be attached to the 
expression V*gr.) We now have 

^\ = VV/>QdS = - V/J/VuQds [by § 9] 

equation (9) § 6 being applied and the centre of the sphere being 
not considered as a singular point since the condition of § 7 is 
satisfied, viz. that Lt r(p_p')=o T^{p — p)uU {p — p) = 0. Now put- 
ting P above = VQ we see that the first expression, viz. VJJJuVQd% 
can be neglected and the second gives 

where Q is the mean value of Q over the surface of the sphere and 
therefore in the limit = Q, Thus equation (19) has been proved. 

When Q has a simple scalar value all the above propositions, 
and indeed processes, become well-known ones in the theory of 
gravitational potential. 

We do not propose to go further into the theory of Potentials 
as the work would not have so direct a bearing on what follows as 
these few considerations. 



SECTION III. 
* Elastic Solids, 

Brief recapitulaticyti of previous work in this branch. 

11. As far as I am aware the only author who has applied 
Quaternions to Elasticity is Prof. Tait. In the chapter on KLine- 
matics of his treatise on Quaternions, ^ 360 — 371, he has considered 
the mathematics of strain with some elaboration and again in the 
chapter on Physical Applications, §§ 487 — 491, he has done the 
same with reference to stress and also its expression in terms of 
the displacement at every point of an elastic body. 

In the former he has very successfully considered various 
useful decompositions of strain into pure and rotational parts and 
so far as strain alone is considered, i.e. without reference to what 
stress brings it about he has left little or nothing to be done. In 
the latter he has worked out the expressions for stress by means 
of certain vector functions at each point, which express the elastic 
properties of the body at that point. 

But as far as I can see his method will not easily adapt itself 
to the solution of problems which have already been considered 
by other methods, or prepare the way for the solution of fresh 
problems. To put Quaternions in this position is our present 
object. I limit myself to the statical aspect of Elasticity, but I 
believe that Quaternions can be as readily, or nearly so, applied to 
the Kinetics of the subject. 

For the sake of completeness I shall repeat in my own no- 
tation a small part of the work that Tait has given. 

* [Note added, 1892. It would be better to head this section " Elastic bodies " 
since except when the strains are assumed small the equations are equally true of 
solids and fluids. I may say here that I have proved in the Proc, R, S, E, 1890 
— 91, pp. 106 et seq.i most of the general propositions of this section somewhat 
more neatly though the processes are essentially the same as here.] 
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Tait shews (§ 370 of his Quaternions) that in any small portion 
of a strained medium the strain is homogeneous and (§ 360) that 
a homogeneous strain function is a linear vector one. He also 
shews (§ 487) that the stress function is a linear vector one and he 
obtains expressions (§§ 487 — 8) for the force per unit volume due 
to the stress, in terms of the space- variation of the stress. 



Strain, Stress-force, Stress-couple. 

12. This last however I give in my own notation. His 
expression in § 370 for the strain function I shall throughout 
denote by x- Thus 

J^Ci)= ft) — ScoV .7] (1), 

where 97 is the displacement that gives rise to the strain. 

Let X consist of a pure strain yjr followed by a rotation q{ )q~^ 
as explained in Tait's Quaternions, § 365 where he obtains both q 
and yjr in terms of x- Thus 

'Xco = qyjrcoq-' (2). 

When the strain is small yfr takes the convenient form ^ where 
'X^ stands for the pure part of x so that 

Xco = ft) — ^Sa)V . 7) — ^^^Soy^i (3)> 

by equation (3) § 2 above. Similarly q{ )q~^ becomes Vd{ ) where 
2^ = Wr], The truth of these statements is seen by putting in 
equation (2) for 3, 1 + ^/2 and therefore for q~^, 1 — ^/2 for y^, x 
and then neglecting all small quantities of an order higher than 
the first. 

13. Next let us find the force and couple per unit volume due 
to a stress which varies from point to point. Let the stress 
function be </>. Then the force on any part of the body, due to 
stress, is 

by equation (9) § 6. Thus tlw force per unit volume = <f>A, for 
the volume considered in the equation may be taken as the 
element d%. 

Again the moment round any arbitrary origin is 
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by equation (9) § 6. The second term on the right is that due to 
the force ^A just considered, and the first shews that in addition 
to this there is a couple per unit volume = V^(l>^ = twice the 
" rotation " vector of <^. Let then ^ be the pure part and Fe ( ) 
the rotatory part of <f). Thus 

Force per unit volume 

= ^A = $A- We (4), 

Couple per unit volume 

= F?<^?=2e (5). 

These results are of course equivalent to those obtained by 
Tait, Quaternions, §§ 487 — 8. 

The meanings just given to r], ;^, p^, yfrj q, <}>, <f> and 6 will be 
retained throughout this Section. In all cases of small strain as 
we have seen we may use '^ or ^ indifferently and whenever we 
wish to indicate that we are considering the physical phenomenon 
of pure strain we shall use yjr, % being regarded merely as a 
function of j^. We shall soon introduce a function ^gf which will 
stand towards (f) somewhat as yjr towards p^ and such that when 
the strain is small 'ut = <f). 

It is to be observed that (jxo is the force exerted on a vector 
area, which when strained is ©, not the stress on an area which 
before strain is eo. Similarly in equations (4) and (5) the in- 
dependent variable of differentiation is p-{-rj so that strictly 
speaking in (4) we should put ^p+,A — Vp^ri^e. In the case of 
small strain these distinctions need not be made. 



Stress in terms of strain, 

14. To express stress in terms of strain we assume any 
displacement and consequent strain at every point of the body 
and then give to every point a small additional displacement St; 
and find in terms of y^ and <^ the increment JJJSwd%Q in the 
potential energy of the body, wd%Q being the potential energy of 
any element of the body whose volume before strain was ds^. Thus 

JJJSwd%Q = (work done by stresses on surface of portion considered), 

— (work done by stresses throughout volume of same 
portion). 
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Thus, observing that by § 12 the rotation due to the small dis- 
placement St) is Vp+j^Sr]/2, we have 

The first of the terms on the right is the work done on the surface 
of the portion of the body considered ; the second is — (work done 
by stress-forces </>p+,,A) ; and the third is — (work done by stress- 
couples 26). Thus converting the surface integral into a volume 
integral by equation (9) § 6 above 

Limiting the portion of the body considered to the element d% we 
get 

Bw = -mSBrj^(t)f,+^V^ (6) 

where m is put for ds/ds^ and therefore may be put by §§ 3a and 3 
above in the various forms 

= s^,^^Sy^^x%xX m 

= 5?.r,C, SKKV^fa (6c) 

= svy;7^S(p + v)Ap+vUp + v)s (6d). 

It is to be observed that since the rotation-vector e of ^ does not 
occur in equation (6), e and therefore the stress-couple are quite arbi- 
trary so far as the strain and potential energy are concerned. I do 
not know whether this has been pointed out before. Of course 
other data in the problem give the stress-couple. In fact it can be 
easily shewn that in all cases whether there be equilibrium or not 
the external couple per unit volume balances the stress-couple. 
[Otherwise the angular acceleration of the element would be 
infinite.] Thus if M be the external couple per unit volume of 
the unstrained solid we have always 

M-|-2m€ = (7). 

In the particular case of equilibrium F being the external force 
per unit volume of the unstrained solid we have 

F + m(^,+,A-Fp+,V€) = (8). 

The mathematical problem is then the same as if for F/m we 
substituted T/iu + Vp+r!^ (M/2m) and for M, zero. In the case of 



§ 14.] ELASTIC SOLIDS. 29 

small strains m may be put = 1. In this case then the mathe- 
matical problem is the same as if for P we substituted P + WM/2 
and for M zero. 

Eeturning to equation (6) observe that 

[which is established just as is the equation x^ = — SooV . (p + v)]> 
so that changing od which is any vector into x^^cd 

Therefore by equation (6) of this section and equation (4) § 3 above 
we have 

We must express the differential on the right of this equation in 
terms of hy^ and hq, the latter however disappearing as we should 
expect. Now ^cd = qyfrcoq'^ [equation (2) § 12] so that remember- 
ing that S . g"' = — q'^Sqq'^ 

= 2VVSqq''^ . %« + qhy^wq'^y 
Sw/m = - 2fif . VSqq-' . f^? - SqSfx''^q~'K 
or since V^(b^= 0, (f> being self-conjugate 

This can be put into a more convenient shape for our present 
purpose. First put for x~^ its value "^'^q'^ ()q and then apply 
equation (6a) § 3 above, putting for the <^ of that equation q~^ ()q 
and therefore for the (f>\ q() q"^. Thus 

Sw = - mSSyfrylr-'^nr^ (9), 

where 'urco = q~^<}) (qcDq~^)q (10). 

The physical meaning of this last equation can easily be shewn. 
Suppose when there is no rotation that (f) = tsr'. Then it is 
natural* to assume — in fact it seems almost axiomatic — that the 
superimposed rotation q()q~^ should merely so to speak rotate 
the stress along with it. Thus if co is some vector area before the 
rotation which becomes ©' by means of the rotation 

q^GT^o) . q"^ = (f>o)\ 

* Observe that we do not make this assumption. We really shew that it is 
true. 
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But <o = gG)g"\ SO that 

Thus we see that nr is as it were </> with the rotation undone. 

Before proceeding further with the calculation let us see what 
we have assumed and what equation (9) teaches us. The one thing 
we have assumed is that the potential energy of the body can be 
taken as the sum of the potential energies of its elements, in other 
words that no part of the potential energy depends conjointly on 
the strains at P and Q where P and Q are points separated by a 
finite distance. This we must take as an axiom. By it we are 
led to the expression for Bw in equation (9). This only involves 
the variation of the pure strain 'yjr but not the space differential 
coefficients of ylr. This is not an obvious result as far as I can 
see but it is I believe always assumed without proof. 

We may now regard w as a function of -i/r only. Therefore by 
equation (7) § 5 above 

therefore by equation (9) of this section 

In equation (6a) § 3 above putting (f) = yjr'^ the right-hand member of 
this equation becomes wfifS^lrfcr^/r'^f. Now putting in equation (6a) 
§ 3 <^ = 'cj^|r"* this member becomes mS^Syjrylr~^tff^ or m/S8'\|rf^"''crf. 
Thus we have 

Now 'DT^|r"* 4- y^'^-rff is self-conjugate. Hence by § 4 above 

m (tsr^lr"^ + ^Ir'^tsr) = 2^QW (11). 

This equation can be looked upon as giving tsr in terms of the 
strain. We can obtain 'cj however explicitly for y^'^'ur is the conju- 
gate of -cj^lr"*. Hence from equation (11) (because ('»/r~''DT -h myjr~^)/2 
is the pure part of CT\|r~*) 

where ^ is a vector to be found. Changing cd into '^(o 

m^co = ^(JWyjroi} -f Vffylrco, 

* By putting d\l/ = u)S{ )w' + w'<S( )a; in this equation, equation (11) can be 
deduced but as this method has ahready been applied in § 4 I give the one in 
the text to shew the variety of Quaternion methods. [Note added, 1892. If we 
use the theorem in the foot-note of § 4, equation (11) follows at once.] 
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Now «r being self-conjugate V^'ut^=^ 0. Hence 

= es^'f^ + ire, 



whence 






...(13). 



or 

by equation (6a) § 3 above. Thus finally 

mvro) = ^(jwylro) + V0ylro) 
or m-BTCD = ^QWyfreo — F^lr© (^ + S^^yjr^J'^Vyjr^^dw^ 

This completely solves the problem of expressing stress in 
terms of strain in the most general case. 

15. 'BTCD 4- Fe'ft), where, as we saw in last section, e is perfectly 
arbitrary so far as the strain is concerned, is the force on the 
strained area cd due to the pure strain ^fr. And again 

(fxo + Ve'o) or q (urq'^eoq) q"^ + Ve'co 

is that due to the strain qylr ( ) q~^ or 'x^. 

To find the force on an area which before strain was (o^ let its 
strained value after yjr has taken place be cd. Then by equation 
(4), § 145 of Tait's Quaternions, mco^ = '>^r(o. Hence 

Required force = 'cjcd + Veto = .^(JWcd^ + FS^o + 'niVe'^^r'^o}^ 

by equation (13) of last section. If the rotation now take place 
this force rotates with it so that the force on the area which was 
originally (o^ is after the strain ;j^ or gi^ {)(f^ 

T^o = 9<Aa^«o7"' + 5'F^«o • ?"' + Vieq'y^-^o)^ , q'') (14) 

where 6 = mqe'q^ and 

.*. e = ^ the stress-couple per un. vol. of unstrained body. 

This force then is a linear vector function of cd^, but in general 
even when 6 = it is not self-conjugate. When both e and the 
rotation are zero we see that the rotation vector of t is ^ given by 
equation (12) of last section. 

The stress-force can be shewn as in § 13 to be tA per unit 
volume of the unstrained body. Thus since the corresponding 
stress-couple is 2e the moment exerted by the stresses on any 
portion of the body round an arbitrary origin is 
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But this moment may also be put in the form 

or Iimr^d%, + JfJVip + v) T.V,d8, + ///F^.rV.ds., 

by equation (9) § 6 above. Comparing these results 

in the notation of next section. This equation may also be 
deduced from equation (15 I) below but not so naturally as above. 

The equations of equilibrium. 

16. We require equation (9) § 14 above to prove the state- 
ment that no space-variations of -^ or g are involved in w. It is 
also required to shew that the fact that q also is not involved 
in ti; is a mathematical sequence of the assumption that the 
potential energy of a solid is the sum of the potential energies of 
its elements. Assuming these facts however we can arrive at the 
equation of stress (11) § 14 in a different way from the above. 
We shall also obtain quite different expressions for ^gf, t, &c. 
and most important of all we shall obtain the equations of 
equilibrium by obtaining t explicitly in terms of the displacement 
and its space derivatives. From § 12 above we have 

Xco = qylrQ)q'\ x^ = '^ (9''^9)> 

XX=r='^ (15), 

(say) as in Tait's Quaternions, § 365. Thus ^cd = V^ScoV^Sp^'p^' 
where p is put as it will be throughout this section for p +rj, the 
vector coordinate after displacement of the point p. From this as 
we have seen in the Introduction we deduce that the coordinates 
of "^ are the A, B, G, a, b, c of Thomson and Tait*s Nat Phil. 
App. C. 

We may as do those authors regard w as a function of ^. 
Thus :— 

= - ssxx^s,aw^- Sx^x^^a^^^' 

By equation (4a) § 3 above, each of the terms in this last 

expression 

= - SBtj^x^q^wV^, 

Bw = - 2SSrj^X*(jwV ^ (15a), 
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Comparing this with equation (6) § 14 and putting in both 
Sr)^ = mSoop where ©', o) are arbitrary constant vectors, we get 

Hence, since a/ is quite arbitrary 

m^ = 2%*a«^x' (1^^)- 

From equation (10) § 14 above which defines 'ur we see that 

mtiT = 2^}r^(JWylr (15c), 

which we should expect since we have already seen that isr is the 
value of (j) when there is no rotation and therefore 'x,=^X='^' 
Now since 

we deduce by any one of the processes already exemplified that 

where of course the differentiations of ^q must not refer to yfr. 
We see then from equation (15c) that 

which is equation (11) § 14. 

Our present purpose however is to find the equations of 
equilibrium. Let co^ be the vector area which by the strain % 
becomes o). Thus* as in last section 

mcoQ = xc^ (15e). 

Further let 2€ be the stress-couple per unit volume of the 
unstrained solid so that 2€/m is the same of the strained solid. 
As we know, e is quite independent of the strain. By § 13 
we see that the force tcDq on the area which before strain was (o^ is 
0o + VeoDJm, Therefore 

To>o = 2x^aw(o^ + ^^X'"'^a (IV)- 

We saw in last section that the force per unit volume of the un- 
strained solid is tA and the couple 2e. Hence 

P + tA = (15^), 

M-h2e = (15/0, 

* See §83 below. 
M. 3 
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are the equations of equilibrium, where P and M are the external 
force and couple per unit volume of the unstrained solid. All that 
remains to be done then is to express t in terms of e and the 
displacement. Putting 

P + V='P (150> 

as already mentioned, we have 

X<o /,/5a,V. (15 j). 

so that by equation (6A) § 3a above we have 

X'-»6, = -3Fp>/^fa,V,V./^.V,V,<Sp>X (Uk). 

Therefore by equation (15/) 

Tfi) 2p/<SV,*a«'a)-3F€F/>,>;<Sfa)V.V,/<SV,V,V,Sp,>,>.'...(150. 

It is unnecessary to write down what equation (i5g) becomes 
when we substitute for t, changing e into — M/2 by equation 
(15A). In the important case however when M = 0, the equation 
is quite simple, viz. 

P = 2/);fifV^^a^A (15m). 

Addition ^o § 16, Dec, 1887 (sent in with the Essay). [The 
following considerations occurred just before I was obliged to send 
the essay in, so that though I thouglit them worth giving I had 
not time to incorporate them in the text. 

It is interesting to consider the case of an isotropic body. 
Here w; is a function of the three principal elongations only and 
therefore we may in accordance with §§14 and 15 suppose it a 
function of a, 6, c or in accordance with § 16 of A, B, G where 

a = -S?t?, b^-S^r^l2, c = ^S^r^/3 (A). 

A = - fiff^f, 5 = - /Sf?>P^?/2, = - S^^'^/S (B), 

Let us use x, y, z, X, F, Z for the differential coefficients of w 
with respect to a, b^ c, A, B, C respectively. Thus 

dw = xda -f ydb + zdc = — 8d'>^^ {x + y^ + Z"^^) f 
as can be easily proved by means of equation (6a) § 3 above. But 

dw = — Sdylr^^i,(jw^. 

^(JW = X + yylr -{- zyjr^ (C^), 

by § 4 above. Similarly we have 

^Qw = X+ Y^-\-Z^'' (D), 
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(Notice in passing that to pass to small strains is quite easy 
for ^(jw is linear and homogeneous inyfr — 1 so that ^(jw = a? + yyjr 
where 3/ is a constant and x + y a, multiple of a — 3.)* 

From equation (C) we see that in equation (12) § 14 is zero 
and therefore that from equation (13) 

mtsr = ,i,(JW^ = ooyfr + yyfr^ + Zyjr^ (E), 

Again from equation (14) § 15 

Similarly from equations (156) and (15/) § 16 we may prove by 
equation (D) that 

im^ = Zxx'+F(xxT + ^(xxT m 

and T = 2Zx + 2 Yxxx + ^^XXXXX + V^x" ( ) (^• 

If we wish to neglect all small quantities above a certain order 
the present equations pave the way for suitably treating the 
subject. I do not however propose to consider the problem here 
as I have not considered it sufficiently to do it justice.] 

Variation of Temperature. 

16a. The w which appears in the above sections is the same as 
the w which occurs in Tait s Thermo-dynamics, § 209, and therefore 
all the above work is true whether the solid experience change of 
temperature or not. to will be a function then of the temperature 
as well as y^. To express the complete mathematical problem of 
the physical behaviour of a solid we ought of course instead of the 
above equations of equilibrium, to have corresponding equations of 
motion, viz. equations (15A), (15i) and (instead of (15gr)) 

2)i5' = P + TA (15n), 

where D is the original density of the solid at the point con- 
sidered. Further we ought to put down the equations of con- 
duction of heat and lastly equations (16/) and (16d) below. 

We do not propose to consider the conduction of heat, but 
it will be well to shew how the thermo-dyn amies of the present 
question are treated by Quaternions. 

* [Note added, 1892. In the original essaj' there was a slip here which I have 
corrected. It was caused by assuming that ^ instead of ^ - 1 was small for small 
strains. In the original I said *' where ?/ is a constant and a; is a multiple of «."] 

3—2 
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Let t be the temperature of the element which was originally 
d%Q and EdS^ its intrinsic energy. Let 

where M, JV" are a linear self-conjugate function and a scalar 
respectively, both functions of yjr and t Here SI{d%Q is the heat 
required to be put into the element to raise its temperature by St 
and its pure strain by Syfr. Now when t is constant SH must be a 
perfect differential so that we may put 

M=t^qf, 

where /is some function of -x/r and t Thus 

BH=^tSSyjf^^af^+^^t (16). 

Now we have seen in § 14 that the work done on the element 
during the increment Sslr, divided by d%Q 

Thus by the first law of Thermo-dynamics 

SE = JSH - mSBylryfr-'^^^ 

where J is Joule's mechanical equivalent. Thus 

JN^dEjdt (16a), 

and ^m (tiT>|r'^ + >|r"* -bt) = ^Qw (166), 

where w = E-Jtf (16c). 

To apply the second law we go through exactly the same cycle 
as does Tait in his Thermo-dynamics, § 209, viz. 

(i|r, t) (f + S^, (^ +Bylr,t+ Bt) (ylr, t + St) (^, t). 

We thus get* 

or Jf=-dwldt (16d), 

the arbitrary function of t being neglected as not affecting any 
physical phenomenon. Substituting for w from equation (16c), 

Jtdf/dt = dE/dt = JN (16c) 

by equation (16a). Thus from equation (16) 

SH=tBf (16/). 



« 



[By assuming from the second law that the work done by the element in the 
cycle, i.e. the sum of the works done by it during the first and third steps is Jdt/t 
multiplied by the heat absorbed by the element in the third step. Note added, 1893.] 
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SO that in elastic solids as in gases we have a convenient function 
which is called the " entropy ". Thus the intrinsic energy EdSo, 
the entropy /d%Q and the stress «r have all been determined in 
terms of one function w of y]r and t which function is therefore in 
this general mathematical theory supposed to be known. 

If instead of regarding w (which with the generalised meaning 
it now bears may still conveniently be called the "potential energy" 
per unit volume) as the fimdamental function of the substance we 
regard the intrinsic energy or the entropy as such it will be seen 
that one other function of yjr must also be known. For suppose / 
the entropy be regarded as known. Then since dw/dt =^ — Jf 

w=W-JJfdt (16g), 

where the integral is any particular one and W is a function of yfr 
only, supposed known. Again 

E = w + Jtf 
or E=W+J(tf-Jfdt) (16h), 

Thus all the functions are given in terms of/ and W. Similarly 
if j& be taken as the fundamental function 

w = t{W'-'fEdtlf) (16i), 

Jf=Elt-\-JEdt/t'-W (16j), 

where as before the integral is some particular one and W is a 
function of yjr only. 

Small strains. 

17. We now make the usual assumption that the strains are 
so small that their coordinates can be neglected in comparison 
with ordinary quantities such as the coefficients of the linear 
vector function ^(jw. We can deduce this case from the above 
more general results. 

To the order considered q = l so that by equation (10) § 14 above 

<^ = -ST. 

We shall use the symbol tsr rather than 4> f^r the same reasons 
explained in § 13 above as induce us to use ^ rather than x- 

Remembering that eo^ and co are now the same and that ^/r may 
be put = l so that V^^(jwylr^= V^^(jw^ = and therefore of 
equation (12) = 0, we have from equation (13), 

OTO) = ^(JWO) (17). 
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Of course we do not require to go through the somewhat com- 
plicated process of §§ 14, 3 5 to arrive at this result. In fact in 
equation (6) § 14, we may put m = 1 and p+,V =V so that 

and therefore by § 3 above 
But by § 5 above 

and therefore by § 4 we get equation (17). 

18. It is convenient here to slightly change the notation. 
For yJTy X we shall now substitute ^/r + 1, j^ + 1 respectively. This 
leads to no confusion as will be seen. 

With this notation the strain being small the stress is linear 
in yjr i.e. ^(jw is linear and therefore w quadratic. Now for any 
such quadratic function 

w = ^Sylr^^aw^/2 (18), 

for we have by § 5, 

dw = — Sdylr^^Qw^. 

Put now yjr = ?i^'. Then because ^Qti; is linear in y^ 

where ,i,(j'w is put for the value of ^(jw when the coordinates of 
yjr' are substituted for those of -x/r. Thus keeping yfr' constant and 
varying ?i, 

whence integrating from n = to n = l and changing yjr' into yjr 
we get equation (18). 

From equation (18) we see that for small strains we have 

W = - Sylr^'ST^/2 (19). 

Now w is quadratic in yfr and therefore also quadratic in tsr, so 
that regarding w as a function of tsr we have as in equation (18) 

w = - S'ST^ ^(lw^/2, 

so that by equation (19) and § 4 above 

t = ^a^ (20). 

All these results for small strains are well-known in their Car- 
tesian form, but it cannot be bias that makes these quaternion 
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proofs appear so much more natural and therefore more simple 
and beautiful than the ordinary ones. 

19. Let us now consider (as in § 16 is really done) w os a, 
function of the displacement. Now w is quadratic in ^/r, and ^Ir is 
linear and symmetrical in V^ and rj^. In fact from equation (3) 
§ 12 above, remembering that the x of that equation is our present 
-x/r + 1, we have 

Therefore we may put 

w==w(vi,V^,rf^,V,) (21), 

where w (a, j3, 7, S) is linear in each of its constituents, is sym- 
metrical in a and ^, and again in 7 and S, and is also such that 
the pair a, )8 and the pair 7, S can be interchanged. [This last 
statement can be made true if not so at first, by substituting for 
w (a, 13, 7, S), w (a, )8, 7, 8)/2 -f w (7, S, a, ^)/2, as this does not affect 
equation (21).] Such a function can be proved to involve 21 inde- 
pendent scalars, which is the number also required to determine 
an arbitrary quadratic function of yjr, since yfr involves six scalars. 

Thus we have the two following expressions for Sw, which we 
equate 

- SB^Ir^^^ = W {Brj^, V„ rj,, VJ + w (^„ V,, S17,, V^) 

or* by § 3 above. 

Now let us put Brj = oo'Scop where a and co are arbitrary constant 
vectors. We thus get 

/SfcoVft) = - 2^ (©', ft), i7i, Vj) = 2Sco'^w (f, ft), i7i, Vj). 

Whence since ft)' is quite arbitrary, 

^ft) = 2?i(;(?,ft),7;,,VJ (22). 

The statical problem can now be easily expressed. As we saw 
in § 14, equations (7) and (8), it is simply 

F + FVM/2 -f tsrA = (22a), 

* [Note added, 1892. Better thus:— by § 3 above, 

- S^^ptm^=2^v {^^pt. f, 7/1, Vi)= - 2S^^^^ptw (fi, r, Viy Vi) 
therefore by § 4, mu}=2l^io (i", w, t/j, Vi) 

for ^w (f, w, 77i, Vi) regarded as a function of w is clearly self conjugate. 
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throughout the mass; and at the surface 

P,-Fra2M/2-i!rC7d2 = (226), 

where P, M are the given external force and couple per unit 
volume and F, is the given external surface traction per unit 
surface. Substituting for -cr from equation (22) 

P + 7VM/2 + 2^w (?, A, ^„ V,) = ) 

p,- vud'ZM/2-2^w(^, crd2,i7„v^)=oj ^ ^* 



Isotropic Bodies. 

20. The simplest way to treat these bodies is to consider the 
(linear) relations between 'gt and '»/r. 

In the first place notice that yfr can always be decomposed into 
three real elongations (contractions being of course considered as 
negative elongations). Thus i being the unit vector in the direc- 
tion of such an elongation, 

^frco = — SeiSioo, 

The elongation — eiSico will cause a stress symmetrical about the 
vector i, i.e. a tension Ae in the direction of i and a pressure Be in 
all directions at right angles ; A and B being constants (on account 
of the linear relation between tsr and sir) independent of the direc- 
tion of i (on account of the isotropy of the solid). This stress 
may otherwise be described as a tension (A -{- B)e in the direction 
of i and a hydrostatic pressure Be, Thus 

-DTft) = — {A -{- B) XeiSico — BcoZe 
= (il + i?) -fo) + BwS^y^i;. 

To obtain the values of A and B in terms of Thomson and 
Tait's coefficients k and n of cubical expansion and rigidity respec- 
tively ; first put 

y^o) = e(o and -cjo) = Skeo), • 
and then put 

ylrcj = VXcofi and 'btg) = 2nV\(t}fi, 

\ and fi being any two vectors perpendicular to each other. We 

thus get 

A+B = 2n, B=:'-{k- 2)1/2), 

whence 

^(o = 2nyfrco - {k - 2n/S) (oS^ylr^ (24). 
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From this we have 

yfra) = m(o/2n + ©/Sff^f {k - 2n/3)/2w, 
but from the same equation 

Equation (24) gives stress in terms of strain and (25) the converse. 

21. We can now give the various useful forms of w for 
isotropic bodies for from equation (19) § 18, 

w = - Sf^7ST};i2. 

Therefore from equations (24) and (25) respectively 

■w n (f^y + i(k- 2n/3) 5??.^?.^??,^ ?, (26). 

- = -ii«)'-4(^-^)«r.-«< (27). - 

Therefore again from § 3 above and from equation (26), 

w = - nS»7,i^V, + ^(k- 2n/3) (fiV,,)' 
or 

+ (k-2n/S)Ss7,v,SV,V, (28). 

Hence from equation (22) 

OTft) = — nSooV , Tj — rN^Srj^o) — (m — n) (oS>7rf (29), 

where m is put for k + n/S. This last could have been deduced at 
once from equation (24) by substituting for -i/ro). 

Thus the equations (23) for the statical problem are 

F+V'7M/2 = nV'7j + 7nVS'7rf (30), 

- P, + VUdXM/2 = nSUdVJ . r, + nV^Srj.Udl 

+ {m-n)UdXSVrf (31). 

We now proceed to apply these results for small strains in 
isotropic bodies to particular cases. These particular cases have 
all been worked out by the aid of Cartesian Geometry and they 
are given to illustrate the truth of the assertion made in the 
Introduction that the consideration of general problems is made 
simpler by the use of Quaternions instead of the ordinary methods. 
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Particular integral of equation (30)*. 

22. Since from equation (30) (P being put for simplicity in- 
stead of P + FVM/2) we have 

nV*77 = P - mVSVri 

we obtain as a particular case by equations (18) and (19) § 10, 

477111; = ///t^ (P - mVS7rj) d%, 

where u has the meaning explained in § 9, and the volume inte- 
gral extends over any portion (say the whole) of the body we may 
choose to consider. To express JJJuVS^T}d% as a function of P put 
in this term u=^—^Up^V^ where p is taken for the p — p' of § 9, 
and apply equation (9) § 6. Thus 

4i7rn7) =fJfuFd% + ^mJJJUpyyS^vd% 

= JJfuFd% - ^mJJJUpV'SVr)d% + a surf. int. 

= JJJuFdi - ^ J^ ffJUpS^Td% + the surf. int. 

for by equation (30) S^T = (m + n)V^S^7f, Now (in order to 
get rid of any infinite terms due to any discontinuity in P) apply 
equation (9) § 6 to the second volume integral. Thus 

^TTurj = ffJuTd% + ^. r ////SfP7 . Upd% + a surf int. 

The surface integral may be neglected as we may thus verify. 
Call the volume integral 4i7rnrj\ Thus 

47rnV V = 47rP + „ , ^^ fJfSFV . V^ Upd% 

4iirmVS^rf' = - JJJSFV . Vud% + . ."^^ , ///>SPV . VSV Upd$, 

n '''''' 2n{m + 7i)'''''' ^ 

so that putting VUp = — 2uvfQ get 

n^^'r}' + mVSVrj' = P, 

whence we have as a pairticular solution of equation (30) rj = 7j' oi 

* [Note added, 1892. For a neater quaternion treatment of this problem see 
Phil. Mag. June, 1892, p. 493.] 
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This is generally regarded as a solution of the statical problem 
for an infinite isotropic body. In this case some law of con- 
vergence must apply to P to make these integrals convergent. 
Thomson and Tait {Nat Phil, § 730) say that this law is that Pr 
converges to zero at infinity. This I think can be disproved by 
an example. Put *Pr = r"^\ where X is a constant vector and a a 
positive constant less than unity. Equation (32) then gives for the 
displacement at the origin due to the part of the integral extend- 
ing throughout a sphere whose centre is the origin and radius R 

m + 37i jB^-^ ^ 

3w (m + n) 1 — a 

Putting R= 00 , 7f also becomes oo . The real law of convergence 
does not seem to me to be worth seeking as the practical utility 
of equation (32) is owing to the fact that it is a particular integral. 

The present solution of the problem has only to be compared 
with the one in Thomson and Tait's Nat Phil §§ 730-1 to see the 
immense advantage to be derived from Quaternions. 

It is easy to put our result in the form given by them. We 
have merely to express /SfPV . Up in terms of P and r^SFV . Vu 
where r is put for the reciprocal of u. Noting that 

Vu = -u^Pj Up=:up, 8¥V ,p = -T, 

we have at once 

8FV .Up=-uF- u'pSTp 

.SfPV . V^ = u^T + ^u'pSFp 

therefore eliminating pSTp, 

STV .Up = - wP2/3 - r'/SfPV . Vu/S, 

.'. 77 = {247rn(m+72)}"V//dS {2(2m+3ii)^P-mr*'/SPV.V2^}...(33), 
which is the required form. 

23. Calling the particular solution rj^ as before and putting 

the statical problem is reduced to finding rj'' to satisfy 

nVy + rnVS^T}" = 

* [Note added, 1892. This is not legitimate since it makes F= oo for r = 0. The 
reasoning is rectified in the PJiil. Mag. paper just referred to by putting F = from 
r=0 to r = b and Fr = r~^\ from r=b to r — cc .] 
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and the surface equation either 

f] + rl' = given value, 
i. e. 7)" = given value, 

or TffUdX = given surface traction, 

i.e. by equation (29) § 21 above, 
nSUdtV . 7i" + uV^St^;' UdX + (m - n) UdXS^v" = known value. 
This general problem for the spherical shell, the only case 
hitherto solved, I do not propose to work out by Quaternions, as 
the methods adopted are the same as those used by Thomson and 
Tait in the same problem. But though each step of the Cartesian 
proof would be represented in the Quaternion, the saving in 
mental labour which is effected by using the peculiarly happy 
notation of Quaternions can only be appreciated by him who has 
worked the whole problem in both notations. The only remark 
necessaiy to make is that we may just as easily use vector, surface 
or solid, harmonics or indeed quaternion harmonics as ordinary 
scalar harmonics. 

Orthogonal coordinates. 

24. It is usual to find what equations (22a) of § 19 and (3) of 
§12 become when expressed in terms of any orthogonal coordi- 
nates. This can be done much more easily by Quaternions than 
Cartesian Geometry. Compare the following investigation with 
the corresponding one in Ibbetson's Math, Theory of Elasticity, 
Chap. V. 

Let Xy y, z be any orthogonal coordinates, i. e. let x = const., 
y = const., -s; = const., represent three families of surfaces cutting 
everywhere at right angles. Particular cases are of course the 
ordinary Cartesian coordinates, the spherical coordinates r, 0, ^ 
and the cylindrical coordinates r, <^, z. Let D^, By, D, stand for 
dififerentiations per unit length perpendicular to the three coordi- 
nate surfaces and let \, //», v be the unit vectors in the correspond- 
ing directions. Thus 

Thus, using the same system of suffixes for the i)s as was ex- 
plained iu connection with V in § 1, 

./.A = D„<l>,\ + D,^<l>,^ + D,,4>,v (34), 

or 

<^A = D« (</)\) + D^ {4>fi) + D, {<\>v) - <t> (Z),\ + D„fi + A")- • -(35). 
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25. Now to put equation (22a) § 19 into the present coordi- 
nates all that is required is to express «rA in terms of those co- 
ordinates. Let the coordinates of cr be PQRSTU. Thus from 
equation (35) we have 

The first thing then is to find Dx\ Dafi, Dy\ &c. 

Letpg, ^3 be the principal curvatures normal to a? = const, i.e. 
(by a well-known property of orthogonal surfaces) the curvatures 
along the lines of intersection o{ x= const., with z = const., and 
y = const, p^y p^ will be considered positive when *the positive 
value of dx is on the convex side of the corresponding curvatures. 
Similarly for q^q^r^r^. Thus for the coordinates r, 9, <f>\ 

JPa = ?8 = 1/^' ?8 = cot eir, q^ = r^==r^ = 0. 
Again for r, (f), z; p^= 1/r and the rest are each zero. 

With these definitions we see geometrically that 

D^ = -- fjq^ - vr^y Da^ = \5,, D^v = \r^ (36). 

Similarly for DyX, Dyfi, DgX &c. Thus 

or 
^A = \ {DJ* + DyU + D,T+ P(p^+p;)- Qp,-Rp, 

+ T(2r, + r,)+U(q,+ 2q;)}+^{} + p{} (37). 

26. The other chiefly useful thing in transformation of coor- 
dinates in the present subject is the expression for the strain 
function yfr in terms of the coordinates of displacement. Let 
u, V, w be these coordinates. Now by equation (3) § 12, remem- 
bering (§ 18) that -^Ir = ^ — 1 we have 

whence - 2yfr\ = XSxD^t) + fiSxDyi] -f vSxD^rj - Dg^Tj, 

But Dx'n = XDxU + fiDxV + vDxW + uDg^ + vDxIJI' + wD^iV 

= X {Dg^u + vq^ + wr^ + /i ( D^v — uq^ + v {D^w — ur^. 

* [Note added, 1892. This is contrary to the usual convention.] 



I (39). 
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Similarly 

DyT} = X (DyU - vp^) + fi (DyV + wr^ + up^ + v (DyW - vr^) 

. • . 2'^\ = 2X (DscV + vq^ + vjr^) + fjL (DyU + D^v -uq^- vp^ 
But with Thomson and Tait's notation for pure small strain 

^e=DxU + vq^ +wr^ ^ 

f = DyV +wr^ -{-up^ 
g = DgW-\-up^ +vq^ 
a = DyW + DzV — wq^ — vr^ 
b = DzU + i)a.w — ur^ — ?^Pg 

Thus we have efgabc in terms of the displacement and we 
have already in equation (24) § 20, which expresses «r in terms of yjr, 
found the values of PQR8TU in terms of e, &c. Finally the 
expression in equation (37) for -crA gives us the equations of equi- 
librium in terms of P, &c. Thus we have all the materials for 
considering any problem with the coordinates we have chosen. 

All these results can be at once applied to spherical and cylin- 
drical coordinates, but as this has nothing to do with our present 
purpose — the exemplification of Quaternion methods — we leave 
the matter here. 

Let us as an example of particular coordinates to which this 
section forms a suitable introduction consider St Venant's Torsion 
Problem by means of cylindrical coordinates. 



Saint-Venani s Torsion Problem. 

27. In this problem we consider the equilibrium of a cylinder 
with any given cross-section, subjected to end-couples, but to no 
bodily forces and no stress on the curved surface. 

* [Note added, 1892. In the Phil. Mag, June, 1892, p. 488, there is a mistake in 
the equation just preceding equation (31) and there are two mistakes in equation (31). 
In the first of these I (2DtU - vsm - jcs w) should be 1l{D^i - tax - /as ?r). In 

equation (31) all the 2's should be dropped.] 
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We shall take r, (f)y z sis our coordinates, the axis of z being 
parallel to the generating lines of the cylinder. Let X, /i, v be 
the unit vectors in the directions of dr, cZ</>, dz respectively and let 

as before. 

We shall follow Thomson and Tait's lines of proof— i.e. we 
shall first find the eflfect of a simple torsion and then add another 
displacement an4 so try to get rid of stress on the curved surface. 

Holding the section z = fixed let us give the cylinder a small 
torsion of magnitude t, i.e. let us put 

r) = Tzr/j, (40), 

for all points for which rz is small. 

The practical manipulation of such expressions as this is almost 
always facilitated by considering the general value of Q(Vj, rj^) 
where Q is any function linear in each of its constituents. Thus 
in the present case 

Q (Vp Vr) = T [zQ (\, /i) - zQ (fi, \) + rQ (p, fi)}. 

[If Q is symmetrical in its constituents, e.g. in the case of stress 

below this reduces to the simple form Q(Vi, ^j) = TrQ(i;, /i).] 

From this we at once see that rj satisfies the equation of internal 

equilibrium 

nVrf + mVSVr) = 0, 

for putting Q (a, ^) = ayS 

Vt) = T (2zv - r\) = tV (z' - rV2) 

so that both SVtj and V^t; = 0. 

Again the value for Q at once gives us the stress for 

OT6D = — nScoV . 7} — nV^ScoTj^ — {m — n) cdSVt;, 

or OTtt) = — nrr (/jlSoov + vSco/x) (41), 

which is a shearing stress nrr on the interfaces perpendicular to 
fjL and V. 

Putting G> = the unit normal of the curved surface we have for 

the surface traction 

'STco = — nrrvScofi. 

In the figure let the plane of the paper be 2j = 0, the origin, 
P a point on the curved surface and OM the perpendicular from 
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on the tangent at P. Thus - rSayfi = OP cos 0PM = PM, PM 
being reckoned positive or negative according as it is in the 




positive or negative direction of rotation round Oz. Thus we see 
that the surface traction is parallel to Oz and = nrPM. 

Hence in the case of a circular cylinder a torsion round the 
axis satisfies all the conditions of our original problem, but this is 
true in no other case. 

The surface traction at any point on the plane ends necessary 
to produce this strain is rsTv^nrrii by equation (41) so that its 
moment round the origin is nrllr^dAy where dA is an element of 
area and the integral extends over the whole cross-section. 

28. Let us now assume a further displacement 

7)=^wv (42), 

where w is a function of r, <j> only, and let us try to determine w 
so that there is still internal equilibrium and so that the stress on 
the curved surface due to w shall neutralise the surface traction 
already considered. 

In the present case 

Thus S^f) = (since v is perpendicular to Vw) and therefore the 
equation of internal equilibrium gives 

V'w = 0. 

Again ma) = — n8(oV . t) — nV^Scorj^f 

or «rft) = — 7i(vSa)Vw '\-VwSg)v) (43), 

a shea,r = nTVw on the interfaces perpendicular to Vw and p. 
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Thus putting cd for the unit normal to the curved surface, the 
present surface traction will neutralize the former if 

ScoVw = — rrScofju = — tSoop (\r) 
= i rSvayV (r»), 

i. e. dw/dn = d (Tr^/2)ld8, 

where d/dn represents differentiation along the normal outwards 
and d/ds differentiation along the positive direction of the bound- 
ing curve. 

We leave the problem here to the theory of complex variables 
and Fourier s Theorem. 

Observe however that the surface traction at any point on the 
plane end = «ri; = nVw by equation (43), and therefore that the 
total couple = - nffrSfjuVwdA = nfj(dw/d(l>) dA. This leads to the 
usual expression for torsional rigidity. 



Wires, 

29. In the following general treatment of Wires some of the 
processes are merely Thomson and Tait s translated into their 
shorter Quaternion forms; others are quite different. The two 
will be easily distinguished by such as are acquainted with Thom- 
son and Tait*s Nat Phil, 

The one thing to be specially careful about is the notation and 
its exact meaning. This meaning we give at the outset. 

The wires we consider are not necessarily naturally straight 
but we assume some definite straight condition of the wire as the 
" geometrically normal " condition. 

The variable in terms of which we wish to express everything 
is s the distance along the wire from some definite point on it. 

Any element of the wire, since it is only slightly strained, may 
be assumed to have turned as a rigid body from its geometrically 
normal position. This rotation is expressed as usual (Tait's 
QuatemionSy § 354) by the quaternion q ; the axis of q being the 
axis of rotation, and the angle of q, half the angle turned through. 

a> is taken so that the rate of this turning per unit length of 
the wire is qtoq'^ so that ta is the rate of turning per unit of length 
when the whole wire is moved as a rigid body so as to bring the 
M. 4 
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element under consideration back to its geometrically normal 
position. Of course o) is a function of q and its derivative with 
reference to s. This function we shall investigate later. The 
resolved part of qoaq^ parallel to the wire is the vector twist and 
the resolved part perpendicular to the wire is in the direction of 
the binormal and equal to the curvature. In fact cd is the vector 
whose coordinates are the /c, X, t of Thomson and Tait's Nat Phil. 
§ 593. When we are given q or (o for every point we know the 
strain of the wire completely, co^ is defined as the naturally 
normal value of co, i.e. the value of a> when the wire is unstressed. 

As usual we take p as the coordinate vector of any point of the 
wire, p like the rest of the functions being considered as a func- 
tion of 5. We shall denote (after Tait, Quaternions^ Chap, ix.) 
diflferentiations with regard to s by dashes. 

We now come to the dynamical symbols. P and M are the 
force and couple respectively exerted across any normal section of 
the wire on the part of the wire which is on the negative side of 
the section by the part on the positive side. 

Finally let X, L be the external force and couple per unit 
length exerted upon the wire. 

30. When the wire is strained in any way let us impose a 
small additional strain represented by an increment h(o in (o and 
an increment he in the elongation at any point. Then the work 
done on the element ds by the stress-force = — SeSTpds and that 
done by the stress-couple = - SqScoq'^Mds, If (as we assume, 
though the assumption is not justified in some useful applications 
of the general theory of wires) P and M to be of the same order 
of magnitude the former of these expressions can be neglected in 
comparison with the latter for Se is a quantity small compared 
with Sq), Now the work done on the element by the stress = the 
increment of the element's potential energy = hwds where w is 
some function of the strain. Hence 

Sw = — SScoq'^'NLq. 
Thus w is a, function of o) only and 

whence we see that 

q''Mq = Js/w (46). 

Notice that M is Thomson and Tait's f, rj, f and q^^Mq their 
KLM (Nat. Phil. §§ 594, 614). 
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31. Now since the strain is small, q'^HLq is linear in terms of 
the strain and therefore in terms of w. Hence we see that w 
is quadratic in terms of a>. Let us then put 

W = Wg (tt) - 0)^, tt) - ft) J + W^ (ft) - (Oq) 

+ a function of temperature only, 

w, and Wj being linear and homogeneous in each of their con- 
stituents. This is the most general quadratic function of ft). 
Now 

Js7w= fWj (f, ft) - ft) J + fw, (ft) - ft),, ?) + fWj?. 

Putting then ft) = ft)o and „Vw = we get ^w^^=0. Operating 
on the last hy So- { ) where <r is any vector we see that w^ = 0. 
Thus putting 

J7w = (t>(co- ft)o), 

where <f> has the value given by the last equation and is therefore 
self-conjugate we get the two following equations 

yj = -8(a>''to,)<f>(a>-a>,)/2 + w, (47), 

[as can be seen by a comparison of the last three equations 
Wq being the function of the temperature] and 

^-^M5 = <^(ft)-ft),) (48). 

When the natural shape of the wire is straight these become 

w = - /Sft)<^ft)/2 + -m;, (49), 

q''Mq = <f>a) (50), 

and when further the wire is truly uniform <f) and w^ are constant 
along the wire. 

32. Assuming the truth of these restrictions let us conceive a 
rigid body moving about a fixed point which, when placed in a 
certain position which we shall call the normal position, has, if 
then rotating with any vector angular velocity &), a moment 
of momentum = cfxo where <f> has the meaning just given. If the 
rigid body be made to take a finite rotation q{) q~^ and then 
to move with angular velocity qcoq'^ its moment of momentum 
will be q<f>(o(f^. Now let a point move along the wire with unit 
velocity and let the rigid body so move in unison with it that 
when the moving point reaches the point s the rigid body shall 
have made the rotation represented by q (§ 29 above). Thus by 
the definition of w and q its angular velocity at any instant 
is qoaq'^ and its moment of momentum therefore q<^(o(f^ or M. 

4—2 
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Now consider the equilibrium of the wire when no external 
force or couple acts except at its ends. In this case P is constant 
throughout and it is easy to see (what indeed is a particular case 
of equation (53) below) that 

M'+F/9'P = (51). 

Interpreting this equation for our rigid body we get as the law 
which governs its motion 

d (vect. mom. of mom.)/cZ^ = — Fp'P. 

Thus the rigid body will move as if acted upon by a constant force 
P at the end of the unit vector p or — since this vector is fixed in 
the body — as if acted upon by a constant force acting through a 
point fixed in the body. From this kinetic analogue of KirchhoflPs 
the mathematical problem of the shape of such a wire as we 
are now considering, under the given circumstances, is shewn to 
be identical with the general problem of the pendulum of which 
the top is a variety. 

33. We will now give the general equations for any wire 
under any external actions. The comparison of the Quaternion 
treatment of this with the Cartesian as given in Thomson and 
Tait's Nat Phil, § 614 seems to me to be all in favour of the 
former. 

The equations of equilibrium of an element ds are with the 
notation explained in § 29 above 

dT + Xds = 0, 

VdpT + (flW + lads = 0, 

or dividing by ds 

P' + X = (52), 

Fp'P + M' + L = (53). 

Operating on the last equation by Vp{ ) noting that p^= — l and 
putting 8p"F = — T we get 

P = />T+F/(M'+L) (54), 

whence by equations (52) and (53) respectively 

X4d[pT+ V(M' + L)}/d5 = (55), 

fifp'(M' + L) = (56). 
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Now by equation (48) above 

5~'Mgr = ^(tt)-tt),) (48). 

Also by the definition of q 

p' = qM'' (57), 

where \ is some given constant unit vector. Finally as we are 
about to prove 

a> = 2Vq''q' (58). 

It is usual in the Cartesian treatment to leave the problem in 
the form of equations equivalent to the above (55) to (58), 13 
scalar equations for the 13 unknown scalars of /?, T, M, q and w. 
We can however as we shall directly in the Quaternion treatment 
quite easily reduce the general problem to one vector and one 
scalar equation involving the four unknown scalars of T and q 
in terms of which all the other unknowns are explicitly given. 

To prove equation (58)* observe that 

{q + dq) <r (g + dq)"^ = gr (o- + Vooads) g~^ 

where a is any vector. The truth of this is seen by noticing that 
(q + dq) ( )(q-\- dq) '* is the operator that rotates any vector of 
the element s + ds from its geometrically normal position to its 
strained position. But we can also get to this final position 
by first in the geometrically normal wire making the small strain 
cods at the given element and then performing the strain of the 
wire up to the point s. The first process is represented on the 
left of the last equation and the second on the right. Thus 
we get 

q^dqa + ad (q~^) . q = Vcjads, 

or • . • d {q"^) = — q'^dqq'^, 

q'^q'a- — a-q~^q = Foxr, 

whence w = 2 Vq'^q, 

Returning to equations (55) to (58) observe that equations 
(57) and (58) give p' and g> as explicit functions of q. Hence 
by equation (48) 

M = ry^ (2 Fj-y -«.)?■' (59). 

* This could be deduced from Tait's Quaternions, § 356, equation (2). His € 
is our quq~^ and his dots our dashes. 
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which gives M also explicitly. Substituting for M and p' in 
equations (55) and (56) we have 

X + ^ \qM-'T + VqM" (t + ^ [# (2 Vq-'q- - a,.) «jr-])| 

= (60), 

Sq^q' (l + ^ [#(2F?-Y -«.)?"']) =0 (61), 

which are sufficient equations to determine T and q, whereupon 
M is given by equation (59), co by equation (58) and p and 
therefore also p by equation (57). 

Spiral springs can be treated very simply by means of the 
above equations, but we have already devoted sufficient space to 
this subject. 



SECTION IV. 
Electricity and Magnetism. 

Electrostatics — General Problem, 

34. Merely observing that all the theorems in integration 
given in the Preliminary and ivth chapters of MaxwelFs treatise 
on Electricity and Magnetism, Part l., are easy particular cases of 
equations (8) and (9) § 6 above, we pass on to the one application 
of Quaternions that we propose to make in Electrostatics. 

This is to find the most general mechanical results arising from 
Maxwell's theory of Electrostatics, and to see if they can be 
explained by stress in the dielectric. This problem as far as I am 
aware has not hitherto in all its generality been attacked though 
the most important practical cases have been, as we shall see, 
considered by Maxwell, Helmholtz, Korteweg, Lorberg and 
KirchhofiF. 

It is necessary first of all to indicate as clearly as possible what 
I take to be Maxwell's theory of Electricity. 

He assumes* all space to be uniformly filled with a certain 
substance called Electricity. Whatever electrical actions take place 

_  

* Prof. J. J. Thomson in his paper on Electrical Theories ^ B.A. Reports, 1885, 
p. 125, does not credit Maxwell with such a definite and circumscribed theory as 
that described in the text, and he is thereby led to find fault with Maxwell's term 
** Displacement" and points out that there is an assumption made with reference 
to the connection between the true current and this polarisation (displacement). 
He says moreover, " It is rather difficult to see what is meant in Maxwell's Theory 
by the phrase * Quantity of Electricity.' " None of these remarks are called for if 
the view I take of Maxwell's theory be correct, and these grounds alone I consider 
sufficient for taking that view. The paper of Thomson's here mentioned I shall 
frequently have to refer to. [Note added, 1892. In the text I have given much too 
rigid a form to Maxwell's theory. What I have called his theory I ought rather 
to have called his analogy. Still I think the present foot-note is in the main just. 
In my opinion it is no more and no less difficult to see what is meant in Maxwell's 
Theory by ** Quantity of Electricity " than by "displacement" since the two are 
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depend on the continued or past motion of this substance as an 
incompressible fluid. If electricity is brought from a distance by 
any means and placed in a given space there must be a displace- 
ment of the original electricity outwards from that space and the 
quantity of foreign electricity is conveniently measured by the 
surface integral of that displacement. 

Dielectrics are substances in which this displacement tends to 
undo itself, so to speak, i.e. the original electricity tends to go back 
to its primitive position. In conductors, on the other hand, there 
is no property distinguishing any imported electricity from the 
original electricity. 

The rate of variation of displacement, whether in dielectrics or 
conductors, of course constitutes an electric current as it is conve- 
niently called. 

We have next to consider a vector at each point of space 
called the electro-motive force, which depends in some way at 
present undefined on the distribution of the displacement in the 
dielectrics, the distribution of currents whether in dielectrics or 
conductors, and on extra-electrical or semi-electrical action, e.g. 
chemical or mechanical. 

If at any point the electro-motive force be multiplied by a 
scalar the medium at the point remaining (except electrically) 
unchanged, the current in the case of conductors and the displace- 
ment in the case of dielectrics is altered in the same ratio. In 
other words the current or the displacement, as the case may be, 
is a linear vector function of the electro-motive force, and the 
coordinates of the linear vector function* at any point depend 

connected by perfectly definite equations. Of course it is wrong to define " dis- 
placement'* as " displacement of quantity of electricity" and then to define "quan- 
tity of electricity" in terms of "displacement," but Maxwell does not seem to me 
even tacitly to do this. Bather he says — the dielectric is polarised ; this polarisa- 
tion can be represented by a vector D ; electrical quantity can be expressed in 
terms of B; the mathematical connections between electrical quantity and B 
are the same as those between quantity of matter in a space and the displacement 
out of that space made by other matter to make room for the given matter ; we 
will impress this useful analogy firmly on our minds by calling B the displacement. 
But I have expressed my present views on the meaning of Maxwell's theory much 
more fully in PhiL Trans. 1892, p. 685.] 

* This frequently recurring cumbrous mode of description must be tolerated 
unless a single word can be invented for "a linear vector function of a vector." 
Might I suggest the term ••Hamiltonian?" Thus we should say that the dis- 
placement is a Hamiltonian of the electro-motive force, the Hamiltonian at any 
point being a function of the state of the medium. 
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solely on the state of the medium (whether fluid, solid, &c., or 
again strained or not) at that point. 

To complete the theory we have to explain how the part of the 
electro-motive force which is a function of the distribution of dis- 
placement and current depends on this distribution. This expla- 
nation is obtained by making the assumption that the electro- 
motive force bears to electricity defined as above exactly the same 
energy relation as ordinary force does to matter, ie. — 

(wk. dn. on electricity moved during any displacement) 

= (total displacement of elect.) x (resolved part of E.M.F. in the 

direction of displacement). 

In the ivth part of Maxwell's treatise he gives complete inves- 
tigations of the mechanical results flowing from this theory so far 
as it refers to currents, but he has not given the general results in 
the case of Electrostatics. Nor has he shewn satisfactorily, it 
seems to me, that the ordinary laws of Electrostatics flow from his 
theory. It is these investigations we now propose to make. 

35. Our notation will be as far as possible the same as 
Maxwells. Thus for the displacement at any point we use D, and 
for the E.M.F. E. From the connection explained in last section 
between D and E we have 

D = irE/47r (1), 

where K at any point is some linear vector function depending on 
the state of the medium at the point. If the medium change in 
any manner not electrical, e.g. by means of ordinary strain K will 
in general also suffer change. 

Let w be the potential energy per unit volume due to the 
electrical configuration. Thus if a small increment SD be given 
to D at all points, the increment jjjhwd% in jjjwd%, the potential 
energy of the electrical configuration in any space, = work done on 
the electricity in producing the change, 

by the relation stated in § 34 existing between E and D. Thus 
limiting the space to the element d% 

Sw = -SIIBJ} (2). 



58 ELECTRICITY AND MAGNETISM. [§ 35. 



Now suppose D = wD' so that by equation (1) E = nB' where 
D' are corresponding E.M.F. and displacement respectively. Thus 

Sw = - nSnSH'I}'. 

Integrating from n=0 to w = 1, and finally changing D', E' into 
jy, E we get 

w = -SDB/2 * (3). 

From this we get 

Stv = - iSfE8D/2 - /gDSE/2, 

so that by equation (2) SESD = /SDSE or by equation (1) 

SBKSB = /SfSEZE. 

Hence (because E and 5E are quite arbitrary) K is self-conjugate 
and therefore involves only six instead of nine coordinates*. 

In electrostatics the line integral of E round any closed curve 
must be zero, for otherwise making a small conductor coincide 
with the curve we shall be able to maintain a current by § 34, and 
so (by the same section) constantly do work on it (i.e. as a matter 
of fact create heat) without altering the statical configuration. 
Hence E must have a potential, say v. Thus 

E = - Vv (4). 

Since in an electrostatic field there is no current in a conductor, 
E = throughout any such conductor and therefore v = const. 

36. The charge in any portion of space is defined as the 
amount of foreign electricity within that space. Thus the charge 
in any space is the surface integral of the displacement outwards. 
Thus if there be a charge on the element d2 of a surface in the 
dielectric this charge = SdSaPa + ^SidSftD^ where a, b denote the 
two faces of the element (so that d^a = — dXb) and in accordance 
with § 1 above dSa points away from the region in which the dis- 
placement is D^. Thus (7 being the surface density 

<T = [Sl}UdtU, (5), 

where the notation [ ]a+u is used for [ ]a + [ ]&• Similarly if there 
be finite volume density of foreign electricity, i.e. finite volume 
density of charge in any space, the charge = --//>SDd2 = - JjJ^SVDdS, 
so that if D be the volume density 

i) = -fifVD (6). 

* We see from this that T>=j,^to or Il=j^^w according as w is looked upon as 
a function of E or D. 
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[The reason for having + Sd^a'^a before and - SdSD here is that 
in the former case we were considering a charge outside the region 
where D^ is considered — between the regions a and b in fact — 
whereas in the latter case we are considering the charge inside the 
region where D is considered. The same explanation applies to 
the sign of — SD Ud% for the surface of a conductor given below.] 

In conductors, as we saw in § 34, the displacement has virtually 
no meaning (except when it is changing and so the phenomenon 
of a current takes place) for the foreign electricity and the original 
electricity are not to be distinguished. Not so however with the 
surface of the dielectric in contact with the conductor. We may 
therefore regard the electricity within the body of the conductor 
as the original electricity so that the charge is entirely at the 
surface. Thus the surface density will be — SD UdX where dX 
points away from the conductor and D is the displacement in the 
dielectric. This may be regarded as a particular case of equation 
(5) D being in accordance with what we have just said considered 
as zero in the conductor. 

37. All the volume integrals with which we now have to deal 
may be considered either to refer to the whole of space or only to 
the dielectrics, as the conductors (except at their surfaces) in all 
cases contribute nothing. The boundary of space will be con- 
sidered as a surface at infinity and all surfaces where either D or XS 
is discontinuous. 

Putting W = JJJwd% we have already found one expression for 

Wy viz. 

2Tr=-///>SfDEds. 

We now give another. By equation (4) 

2W = JJJSDVvd% 

= JJvSdXB - JJJvSVBds, 
by equation (9) § 6 above. Thus by equations (5) and (6) § 36, 

2W = JJvads + JJJvDd% (7), 

where ds is put, as it frequently will be, for an element of surface, 
i.e. Td%. The value of W which we shall use* is obtained by 
combining these two, viz. 

W = JJv<Tds-^JfJvDd% + ^fJfSl}Iid% (8). 

* This is for the general case following the example of Helmholtz in the par- 
ticular case when K reduces to a single scalar. See Wiss. Abh. vol. i. equation {2d), 
p. 805. The method adopted in the following investigation is also similar to his. 
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So far we have merely been shewing that all the above results 
of MaxwelFs flow from what in § 34 has been described as his 
theory. We now proceed to the actual problem in hand which is 
proved from these results however they may be obtained. I may 
remark that some such investigation as the above seems to me 
necessary to make the logic of Maxwell's treatise complete. 

38. Suppose now that W is the potential energy of some 
dynamical system extending throughout space. Let us give to 
every point of space a small displacement Siy vanishing at infinity 
and find the consequent increment 8 W in W. If this can be put in 
the form 

8Tf = -///5S,,<^V.ds (9), 

we shall have the following expression for P the force per unit 
volume due to the system 

P = <^A : (10), 

and the following expression for Pg the force per unit surface at 
any surface of discontinuity in <^ 

T, = -\<j>Udt\^„ (11), 

the notation being the same as in equation (5) § 36. 

Moreover if be self-conjugate the forces both throughout the 
volume and at surfaces of discontinuity are producible by the 
stress <f> as can be seen by § 13 above. [Compare all these state- 
ments with § 14 above.] 

For proof, we have by equation (9) § 6 

S F = - JJJ8Br,,<f>^,d» = -JJS^<t>dX + mBv<j>yA> 

where of course the element dS is taken twicey i.e. once for each 
face. But 

S TT = — (work done by the system P, P^ of forces) 

where the element ds is taken only once. Equating the coefficients 
of the arbitrary vector St) for each point of space we get the 
required equations (10) and (11). 

39. ^We must then put BW where W is given by equation (8) 
in the form given in equation (9). 

We must first define S when applied to a function of the posi- 
tion of a point. Suppose by means of the small displacement Br) 
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any point P moves to P\ Then Q being the value at P, before 
the displacement &;, of a function of the position of a point, Q + SQ 
is defined as the value of the function at P' after the displace- 
ment. Thus even in the neighbourhood of a surface of discon- 
tinuity hQ is a small quantity of the same order as S^. 

Now the charge within any space, that is the quantity of 
foreign electricity within that space will not be altered by the 
strain. 

.-. S(i)d»)=0, h{(Tds) = (12). 

To find hV we have 

S {dp + Up) (7 + S7) . = ScJpV. 

or, since hdp = — Sdp^ . Br), 

SdpSV = SdpV^SBrj.V, 

whence BV =V^Shvy (13). 

The part of STT depending on the first two terms of equation 
(8) is by equation (12) 

JJJBvDd% + JJBvads 
= - JJJBv8Vl}d% + JfSvSdll} [by equations (5) and (6) § 36] 
= JJJ8DVBvd% [by equation (9) § 6]. 

Noticing that Bd% = - d^S^Brj and that 47rSD = KBH + BKH we 
see that the last term in equation (8) contributes 

JJJ8DBBd% - i JffSl}B8VBvd% + (Stt)-' JJJSEBKBd%. 
Combining the last result with the first term of this we get 

fffSB {VBv - B (Vi;)) d» = - ffJSl}BV.vd% 

^-•fffSDV^SBvyvd% 

= JffSI}V^SEBv,d%. 
Thus we have 

g TT = - fffSBff, (i V^/SDE - ESV^D) d% + (87r)-^ JffSEBKBd$. . .(14). 

40. Now the increment BK in K is caused by two things 
viz. the mere rotation of the body and the change of shape of the 
body. Let us call these parts BKr and BKg respectively. 

First consider BKr. Suppose the rotation is e so that any 
vector which was co becomes thereby a> + Veco, Thus the result 
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of operating on o) + Few by K + hK^ is the same as first operating 
on 0) by ^ and then rotating. In symbols 

whence SKrCo = VeKco -K( Veco). 

Thus SBBKr^ = iSfEe^E - SEK ( FeE) 

= 2>Sf€ZEE = 87r/SfeDE, 
whence giving e its value FVS17/2, 

(Stt)-* SaSS^^ = /SfVSi;Fl>E/2. 

Substituting in equation (14) 

SW= - JJfSSf)^ {FV,DE/2 -ESV^D} d% + (Sirr JJJSllSKsEdS, 
or 

S TT = J ///fif&7, FDV^Eds + (87rr UfSESK^ds (15). 

It only remains to consider SK^* K is Si function of the pure 
strain of the medium and SKg is the increment in K due to the 
increment in pure strain owing to St). Calling this increment of 
pure strain Syfr so that by equation (3) § 12 above 

we have 

Sir, = -ss^?^a2?.^2 (i6)> 

by equation (7) § 5 above. This gives by equation (5) § 3 

BKs = '-8Sv,^a^V^,K, (17), 

SESK^ = -8Sv,^(j2'^,SEK^-B (18). 

Now by equations (1) and (3) § 35, 

Sttw = - SERB, 

so that t^ is a function of the independent variables E, sjr (because 
iT is a function of '^jr). Therefore 

(87r)"'/SfESir^ = /SSt^^^q^V^ (19)- 

This equation might have been deduced at once thus 

(S7r)-'8BBKsB = /SfS>|rf^a«^?= SBv.^'DwV^, 

but equation (17) is itself of importance so the above proof is 
preferable. 

Thus finally from equation (15) we get 

S W = JJJSBv, ( FDV,E/2 + ^q^V,) ds (20). 
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We therefore have for <^ in equations (10) and (11) § 38, 

^ft) = - Vl}(aB/2 - ^(jwo) (21). 

This is a self-conjugate function so that as we saw in § 38 it is 
a stress which serves to explain forces both throughout the volume 
of the dielectric and over any surfaces of discontinuity in D or 



The force in particular cases, 

41. Let us first consider that part — FDAE/2 of the force 
(equations (10) § 38 and (21") § 40) which does not depend on the 
variation of K with the shape of the body. 

Suppose our dielectric is homogeneous and electrically isotropic 
so that j&T is a simple constant scalar. In this case 

47rD = -ZVv (22), 

by equations (1) and (4) § 35. Therefore by equations (5) and 

(6) § 36, 

47rD = irV»v (23), 

47r<r ^-KiSUdtVv^a^i, (24). 

From these we at once get by the theory of potential that 

Kv = JJfuDd% + Jfuads (25). 

From this we know by the theory of potential that at the surface 
where the charge a resides Vv is discontinuous only with regard 
to its normal component and at all other points is continuous. 

Thus 

and by equation (24) oo = ^tira-jK so that 

(47r)-^Z {Vva - Vt;^) = aUdta. 

whence (4>7r)-'KUdta (Vv^ - Vv^) = -0- (26). 

Now the force P per unit volume is 

- FDAE/2 = - ZyVvAV-y/STT 

or T = -DVv (27), 

* As far as I am aware nobody has hitherto attempted to find the electrical 
forces much less the stress except in the case when D is parallel to E i.e. the 
dielectric is electrically isotropic when unstrained. The particular results contained 
in § 45 below have been obtained by Eorteweg, Lorberg and EirchhofF as is stated 
in Prof. J. J. Thomson's paper (p. 155) referred to in § 34. 
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and the force per unit surface P« is by equation (11) § 38, 
i [VBUdlHU, = {87rrK[r^vUd^Vv]a+b 

= (SirrKVi^VaUdta^Va ~ Vvj^Udta^V^) 

whence by equation (26) 

P« = -Jo-[VvW (28). 

Thus we see that Maxwell's theory as given in § 34 above 
reduces to the ordinary theory when K is a, single scalar. In fact 
two particles containing charges e e' apparently repel one another 
with a force ee/Kr^ where r is the distance between them, for by 
equations (25), (27) and (28) the force in any charged body is that 
due to a field of potential v given by 

Kv = l.ue (29). 

42. If the medium when strained remain electrically isotropic 
^(jK as well as K must be a simple scalar. Thus with Thomson 
and Tait's notation for strain, which makes the coordinates of 
^, e,f, g, a/2, 6/2, c/2 we have 

^air = dK/de = dZ/d/= dK/dg 
dKjda = dK/db = dKjdc = 0. 

Therefore iT is a function of e +/+ g only, i.e. of the density (m) 
of the medium. Thus because 

de + df+ dg = — dm/m = — d log m 

we get ^(jK = — dK/d log m = — k suppose. Hence 

^d^i; = (SttY'SE (dK/d log m) E = ^EVStt. 

Thus the force — ^Qt^A [equations (21) § 40 and (10) § 38] result- 
ing from the change of K with pure strain is in the case we are 

now considering 

-.A;VE787r (30), 

and is* therefore, according as A; is positive or negative, in the 
direction of or that opposite to that of the most rapid increase of 
the square of the electromotive force. Thus even in the case of a 
fluid dielectric which has no internal charge but which forms part 
of a non-uniform field of (electromotive) force the surfaces of equal 
pressure and therefore the free surface will if originally plane no 
longer remain so. 

* This is the same result as Helmholtz*a on the same assumption Wists. Ahh, 
I. p. 798. 



§ 44.] ELECTRICITY AND MAGNETISM. 65 

Nature of the Stress. 

43. We have seen that the stress which serves to explain the 
electrostatic forces is that given by equation (21) § 40, viz. 

<j)(o = — FDft)E/2 - ^(jwd) (21). 

Let us first consider the part — 7Da)E/2 which does not 

depend on the variation of K, Putting eo first = Ul> and then 

= J7E we get 

^[7B = !rD!rE. trE/2, 

(l>UB = TBTH . UJ}/2. 

Therefore putting w first = any multiple of ?7D + tTE and then 
= any multiple of Z7D — UH we get 

(fxo = TJ}TBcol2, 

(l>(o = - Tl}TIia}/2, 
Lastly, since 

- FDoE = a)/SDE - DjSfa>E - BSeoJ}, 

we see that if we put eo = any multiple of FDE 

(fxa = awSiDE/2 = -- cdw. 

Thus we see that the stress now considered is a tension along 
one of the bisectors of D and E (the bisector of the positive direc- 
tions or the negative directions of both) = 2'DrE/2, an equal 
pressure along the other bisector and a pressure = w perpendicular 
to both these directions. When D is parallel to E this at once 
reduces to Maxwell's case, viz. a tension in the direction of E and 
a pressure in all directions at right angles each = w. 

44. We have now to consider the other part of the stress, viz. 

or <t>a> = ^(Ji^SBK^E/Stt (31). 

If we assume that if is a function of the density (m) of the 
medium only we shall have 

dK/de = dK/df= dK/dg = - dK/d log m = - k, 

say, and 

dK/da = dK/db = dK/dc = 0, 

as in § 42. Here however k is not in general a mere scalar but a 
self- conjugate linear vector function. We have then in this case 

<^ft) = - coSHkH/STT, 
M. 5 
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which is a hydrostatic pressure or an equal tension in all direc- 
tions according as fifEiE is positive or negative. In this case the 
36 coordinates of 4,(liCi) . K^ reduce to the 6 of i for each point of 
space. 

A more general assumption is that SKg (§ 40) depends only on 
the elongations in the directions of the principal axes of K. 
Taking i,j, k as unit vectors in these directions we again have 

dK/da = dK/db = dK/dc = 0, 

and thus <f>i = iSE (dK/de) E/Stt, 

and similarly for j and k, so that the principal axes of the stress 
now considered are the principal axes of K. 

45.* The most natural simple assumption for solid dielectrics 
seems to me to be that the medium is electrically isotropic before 
strain, and also isotropic with regard to the strain in the sense 
that if the strain be, so to speak, merely rotated, SKg will suffer 
exactly the same rotation. We may treat this problem exactly as 
we did (§ 20) that of stress in terms of strain for an isotropic solid. 
Thus splitting up Byjr into its principal elongations, i.e. putting 

Byjrco = — XBeiSio), 
we shall get, as in § 20, 

SKgO} = — (a — /3) XBeiSico + ^(o%he 

But SKgCo = ■- SSslr^^Q^^ , K^co by equation (16) § 40, so that from 
equation (31) 

- STT^SfSilrfc^f = (a - /3) >SESi|rE - ^H'S^Sf^, 

whence we see by § 4 above that 

^0) = {(a - y8) E^SEo) + /3EX/87r (32), 

which consists of a pressure in the directions of the lines of force 
= — aE^/STT and another pressure in all directions at right angles 

= - ^EVStt. 

* For references to former proofs of this see foot-note to § 40 above. 
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Magnetism. 

Magnetic potential, force, induction. 

46. We now go on to the ordinary theory of magnetism ; and 
here we shall merely follow Maxwell in his General Theory, so as 
to give an opportunity of comparing Quaternion proofs with Car- 
tesian, as we have already done in Elasticity. 

We shall not consider in detail the effect of one small magnet 
upon another, as this has already been done by Tait. In connec- 
tion with this I am content to remark that I think the treatment 
of this problem can be made somewhat simpler than Tait's by 
means of potential. 

Suppose we have a pole — m at and a pole + m at 0' where 

00' is small. Calling the vector from to 0' 00\ let us call the 

vector mOO' fi, so that fi is the vector magnetic moment of the 
magnet. The potential of — m at any point P is — mUy where u as 
usual = PO"^. Similarly the potential of + m is mu\ where 
u' = PO^'^. Therefore the potential of the magnet 

= m(u — u) = mSO O^u, 

where of course P is the variable point implied by V. Thus the 
potential of a small magnet //, at any point = S/jl^u. 

Hence the potential of any magnet whose magnetic moment 
per unit volume at any point is I is 

n = S'^fJfuld% = -JfJSlVud% (33), 

according to the convention of § 9 above. By equation (9) § 6 this 

may be put 

n = - fJnSldt + JJJuSVldS (34), 

which shews that we may consider it due to a volume density 
fifVl and a surface density — 81 UdX * of magnetic matter, the sur- 
face density occurring wherever there is discontinuity in I. 

By again considering the poles m and — m of the small magnet 
/Lt we see that its potential energy when placed in a field of mag- 
netic potential fl is — SfjuVn, whence just as we obtained equation 

* [Note added, 1892. More generally and better -[SlUdl.]^^^.] 

5—2 
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(33) we now see that the potential energy (W) of any magnet in 

such a field is 

W = -fffSlVnd% (35), 

or F=-J/ft>SfIrf2 + ///n>SVld» (36), 

by equation (9) § 6 above, so that the potential energy is just the 
same as it would be for the imaginary distribution of magnetic 
matter. 

47. The force (H) on a unit magnetic pole at any point 
external to the magnet is given by 

H = - Vft = - VS^fJfuld% = - V'JfJuld% + V rVfffuld% 

= - 4^1 + VJJJVlVud% = - 47rl + VA, 
where 

A=fffVlVvd% (37). 

Thus we see that for all external points H = VA, so that A is 
called the vector magnetic potential. [It is to be observed that 
since VA = H + 47rl = a vector, >Si VA = 0.] V A is called the mag- 
netic induction and for it we use the single symbol B so that 

B = VA (38). 

Thus iSfVB = iSfV^A = (39), 

and also by the equation for H just given 

B = H + 47rI (40). 

This is not the way in which Maxwell defines the magnetic 
force and induction, but he shews quite simply (Elect, and Mag. 
^ 398 — 9) that his definition and the present one are identical. 
This can be shewn as easily without analysis at all. 

48. Where I is discontinuous both H and B are also discon- 
tinuous. From the surface density view we gave in equation (34) 
we see that, just as we have the expression, given in § 41 for 
Vva—^Vby so now 

llb-l^a=^-^'^Udla[SlUdt]a+b (41), 

SO that the discontinuity in H is entirely normal to the surface of 
discontinuity. Further from this equation we have 

«f (Hft - H^) dXa = 47r [8ld%]a+b, 
i.e. [S (H + 47rl) dX]a+b = 0, 

or [SBdXUb = (42), 

so that the discontinuity in B is entirely tangential. 
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From this equation we see that for any closed surface whatever 
whether it include surfaces of discontinuity in I or not 

JJSdtB = 0. 

For adding these surfaces to the boundary of the inclosed space, 
in accordance with § 7 above, we see by equation (42) that they 
contribute zero to the surface integral ; but the total surface inte- 
gral is by equation (9) § 6 JJJSVBd% = by equation (39). 



Magnetic Solenoids and Shells, 

49. A magnet is said to be solenoidal if the imaginary mag- 
netic matter of equation (34) is entirely on the surface. Thus 
for a solenoidal distribution 

^VI = (43). 

In this case the potential is by equation (34) 

n = - jjuSldt (44). 

50. A simple magnetic shell is defined as a sheet magnetised 
everywhere normally to itself and such that, at any point, the 
magnetic moment per unit surface is a constant called the strength 
of the shell. 

Calling the strength we have for the potential energy at aoy 
point by equation (33) 

n = ^(l>fJSdt^u (45). 

Now — Sdl^Vu is the solid angle subtended by the element dS at 
the point considered, so that 

^ = <f> X (solid angle subtended by shell at point) (46). 

Thus if P be a point on the positive side of the shell and P' a 
point infinitely near P but on the negative side 

Potential at P -- Potential at P' = 47r^, 

or what comes to the same thing 

- / SUdp = 47r(/>. 
J p 

This integral may be taken along any path, e.g. along a path 
which nowhere cuts the shell. The same integral is true if H be 
the magnetic force due to a whole field of which the shell is only 
one of several causes, for the part contributed by the rest of the 
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field is zero on account of the infinite proximity of P and P*. 
For future use in electro-magnetism observe that this statement 
cannot be made if for H in the integral be substituted B. 

51. The condition that any magnet can be divided up into 
such shells is at once seen to be that I can be put in the form 

I = V<^ (47), 

where <f> is some scalar. 

In this case the potential is by equation (33) 

fl = - ffJSlVudS = - fJJ8V<t>Vud%, 

or by equation (9) § 6 

11 = - fJ<t>S(W7u + fJJ<t>V'ud%, 

i.e. a = - JJ<t>8dXVu -{- 4i7r(t> (48). 

Remarking that the solid angle again occurs here it is needless to 
interpret the equation further. By equation (37) we have for the 
vector potential 

A=JJJVlVvd% =JJJW<l>Vmi%, 

or by equation (9) § 6 

A=^JJ(l>rdtVu (49). 

52. The potential eoergy of a magnetic shell of strength ^ 
placed in a field of potential fi is of importance. We see by 
equation (35) that it is 

If then the magnets which cause 12 do not cut the shell anywhere, 
so that — VII = VA, we shall have 

W = (l>JJSdXVA = (l>JJSdtB (50), 

or W= ^JSdpA (51), 

by equation (8) § 6. 

Suppose now that A is caused by another shell of strength <!>'. 
Then by equation (49) 

A = <l>'JJVdt'Vu = <f>'Judp\ 

by equation (8) § 6. Thus finally the potential energy M of these 
two shells is given by 

M = <t>(}>' JJuSdpdp (52). 

53. The general theory of induced magnetism when once the 
proposition (given in equation (42) § 48) that [SdX'B]a+b is zero is 
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established, is much the same whether treated by Quaternion or 
Cartesian notation. We shall therefore not enter into this part of 
the subject. 



Electro-magnetism. 
General Theory. 

54. We now propose to prove the geometrical theorems con- 
nected with Maxwell's general theory of Electro-magnetism by 
means of Quaternions. 

We assume the dynamical results of Chaps. V., VI. and VII., 
and the first six paragraphs of Chap. VIII. of the fourth part of 
his treatise. 

These assumptions amount to the following. Connected with 
any closed curve in an electro-magnetic field there is a function 

p = -fSAdp (53), 

where A is some vector function at every point of the field. The 
function p has the following properties. If any circuit be made 
to coincide with the curve the generalised force acting upon the 
electricity in the circuit is 

E^-p (54). 

Again, if there be a current of electricity 7 flowing round this 
circuit, the generalised force X, corresponding to any coordinate x 
of the position of the circuit due to the field acting upon the con- 
ductor, is 

X = jdpldw (55). 

00. The first thing to be noticed is that p can be transformed 
into a surface integral by equation (8) § 6 above. 

Thus p = -fJSBdl (56), 

where . B = FVA (57), 

so that SVB = 0. 

Next we see by the fundamental connection (§ 34 above) 
between the e.m.f. E and electricity, that E must equal the line 
integral of E round the circuit, or 

E^-'JSEdp (58). 
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We are now in a position to find E in terms of A and B, Le. 

of A. 

56. The rate of variation of p is due to two causes, viz. the 
variation of the field (A) and the motion of the circuit (p). In 
the time it then there will be an increment SA in A and an incre- 
ment MX in d% to be considered. Thus 

pit = "JSiAdp - ///SBScZS. 

[This amounts to assuming that JSiAdp = ffSilBdZ, which of 
course is true by equation (8) § 6.] Now when the circuit changes 
slightly we may suppose the surface over wliich the new integral 
extends to coincide with the original surface and a small strip at 
the boundary traced out by the motion (pit) of the boundary. 
Thus Sd2 is zero everywhere except at the boundary and there it 

= V(pit) dp, 
so that pit = - JSiAdp + it JSpBdp, 

whence dividing by it 

P=J8i-A+VpB)dp (59), 

but by equations (54) and (58) p =fSHdp. Thus 

E = -A+F/i)B--V^ (60), 

where yfr is a. scalar and -- Vyjr is put as the most general vector 
whose line integral round any closed curve is zero. 

57. We now come to the mechanical forces exerted on an 
element through which a current O per unit volume flows. 

We see by equation (55) that the work done by the mechanical 
forces on any circuit through which a current of magnitude 7 flows 
in any small displacement of the circuit equals 7 x the increment in 
p caused by the displacement. Give then to each element dp of 
the circuit an arbitrary small displacement ip and let P' be the 
mechanical force exerted by the field upon the element. Thus as 
in last section 

-JSF'ip = yip = - yfJSBidl, = - yfSipdpB. 

Thus the force P' on the element dp is yVdplB, But we may 
suppose this element to be an element d^ of volume through which 
the current O flows. Thus if for ydp we write Cd%, and for P', 
PdS, where P is the force per unit volume exerted by the field, we 



get 



P= FOB (61). 
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58. So far we Have been able to go by considering the electric 
field as a mechanical system, but to go further (as Maxwell points 
out) and find how B or A depends on the distribution of current 
and displacement in the field we must appeal to experiment. It 
has been shewn by experiment that a small circuit produces 
exactly the same mechanical eflfects on magnets as would a small 
magnet, at the same point as the circuit, placed with its positive 
pole pointing in the direction of the positive normal to the plane 
of the circuit when the positive direction round the circuit is 
taken as that of the current*. Moreover the magnetic moment 
of the magnet which must be placed there is proportional to the 
strength of tbe current x the area of the circuit. Further, the 
effect of this circuit upon other such small circuits is the same as 
the mutual effects of corresponding magnets. We have now only 
to consider a finite circuit as split up in the usual way into a 
number of elementary circuits to see that a finite circuit will act 
upon magnets or upon other circuits exactly like a magnetic shell 
of strength proportional to the strength of the current and 
boundary coinciding with that of the circuit. The unit current 
in tbe electro-magnetic system is so taken as to make this propor- 
tionality an equality. 

The one difference between the circuit and the magnetic, shell 
is that there is no discontinuity in the magnetic potential in going 
round the circuit, so that by § 50 above the line integral of H 
round the circuit will be 47r x the strength of the current. In 
symbols 

fSHdp = 4!7rJfSCdl 
for any curve, so that by equation (8) § 6 

47rO= FVH (62), 

whence SVC = (63), 

which of course is a direct result of our original assumption that 
electricity moves like an incompressible fluid. Maxwell tacitly 
assumes this by making the assumption that only one coordinate 
is required to express the motion of electricity in a circuit. 

* See § 1 above for the convention with respect to the relation between the 
positive side of a surface and the positive direction round its boundary. Hitherto 
we have had no reason for choosing either the right-handed or the left handed 
screw as the type of positive and negative rotation. But to make the statement in 
the text correct we must take the former. 
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59. We are now in a position to identify the B we are now- 
using with the magnetic induction for which we have already used 
the same symbol. 

We see by equation (50) § 52 that the mechanical force on tlie 
shell corresponding to any coordinate x is 

-'<f>dS5SdtB'ldx, 

where B' is the magnetic induction ; and by equations (53) and 
(55) that the force on the corresponding electric circuit is 

- (f>dJJSdXB/dx, 

therefore B = B' wherever there is no magnetism. And where 
there is magnetism B is not = H for £iVB = 0, as we have seen. 
Thus B = B' at all points. In other words the two vectors are 
identical and we are justified in using the same symbol for the 
two. 

This practically ends the general theory of electro-magnetism. 
We content ourselves with one more application of Quaternions 
in this subject. We give it because it exhibits in a striking 
manner the advantages of Quaternion methods. 



Electro-magnetic phenomena explained by Stress, 

60. In § 46 equation (35) we have seen that the potential 
energy of a magnetic element = ld% in a magnetic field is SlH.d% 
when H has a potential. Maxwell assumes that the same ex- 
pression is true whether H have a potential or not. Assuming 
this point* with him we can find the force and couple acting on 
the medium and a stress which will produce that force and 
couple. The force and couple due to the magnetism of an element 



* 



I do not defend the legitimacy of this assumption. It seems to me bold to 
assume that a magnet possesses any such thing as potential energy in a field which 
has no potential. If we assume HE and its derivatives to be continuous throughout 
our typical element d^ of volume containing a great number of molecules (both 
material and magnetic) the force on a magnetic molecule fi consisting of two poles 
is - Sfj.\. HE and the force per unit volume - Sl^. H, which is only identical with the 
expression - ^^HlH.^ obtained below when K has a potential. With this expression 
a stress cannot be found that produces the force. If, however, S and its derivatives 
be not assumed continuous in this manner the force on the magnet /x is quite inde- 
terminate whether the magnetic pole or the molecular current view of magnetism 
be taken, unless it be specified in what way the poles and currents are distributed in 
the element of volume. 
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is obtained by giving the element an arbitrary translation and 
rotation and assuming that the work done by this force and 
couple = the decrement in the potential energy of the element. 
Thus the force per unit volume is — V^^IHi for the decrement in 
the potential energy due to a small translation 8p is SSpV^SlRi. 
Similarly the couple M is given by 

M= 71H (64) 

for the decrement — SMSeo in the potential energy due to a small 
rotation Bco is —SVB(ol,'H, = — SS(oVtH., The total force P per 
unit volume is the sum of that just given and that given by 
equation (61), so that 

P=FOB-ViSIHj. 

Therefore by equation (62) 

47rP = FFV^H . B - 47rV^fifIHj 

= - H,/SViB + V^/SBH, - 4>7rV^Sni^ 

= -HififV^B + VjSHHj. 

Now SVB = so that H,fifV,B = H^AB, and again 

V^SHH, = V(H0/2, 

so that P = <^A (65), 

where Bttc^o) = - 2HSa>B + a>H' (66). 

From this we get 

Stt V^(f>^=2 FBH = Stt FIH, 

so that M=F?05' (67). 

From these two results (equations (65) and (67)) we see by 
§ 13 above that the stress will produce all the mechanical effects 
of the field. 

This stress, as can be seen by giving a> the required values in 
equation (66), is one of pressure — IV /Sir in all directions at right 
angles to B and of tension — /SH (2B — H)/87r in the direction of 
H. When there is no magnetism H = B so that this pressure and 
tension become equal and their directions at right angles to and 
along B respectively. In fact we then have 

87r</)ft) = - Ba)B (68). 



SECTION V. 
Hydrodynamics. 

61. In the applications I am about to make in this I have 
practically nothing new to shew except the utility of Quaternion 
methods in the general theory of Hydrodynamics in all its parts. 

I therefore take a treatise (Greenhiirs article in the Encyc, 
Brit on this subject) and work out the general theory on the lines 
of the treatise. This is more necessary than at first sight it would 
seem, for I believe mathematicians who have studied Quaternions 
are under the impression that the method does not lend itself 
conveniently to the establishment and treatment of such equations 
as the Lagrangian and those of Cauchy's integrals. With our 
meaning of V, however, the Quaternion treatment of these equa- 
tions is as much simpler than the Cartesian as in the case of the 
Eulerian equations. 

Down to equation (13) below the subject has already been 
handled by Prof. Hicks* in his Quaternion treatment of Strains 
and Fluid motion {Quart Journ. Math. xiv. [1877] p. 271). I do 
not hesitate to go over the ground again as my methods are 
dififerent from his. 



Notation, 

62. For the vector velocity at any point we shall use <r, for 
the density D, for the force per unit mass P, for its potential, when 
it has one, v, and for the pressure p. For time-flux which follows 
a particle we shall use djdt or the Newtonian dot, and for that 
which refers to a fixed point of space 9/9^f. 

* [Note added, 1892. Prof. Tail's name ought to be added to Prof. Hicks's.] 
+ [Note added, 1892. I am aware that this is contrary to the usual English 
custom, but that custom — of interchanging the meanings of djdt and dldt as given 
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Thus 

d/dt=^d/dt-8aV (1). 

This equation is given on p. 446 of Greenhill's article already 
mentioned. In future we shall refer to this article simply as 
" Greenbill's article." 



Elder's Equations. 

63. To find the equation of continuity, with Greenhill, we 
merely express symbolically that the rate of increase of the mass of 
the fluid in any space equals the rate at which it is flowing 
through the boundary. Thus M being the mass in any space, 

dMjdt^SJDSadt (2). 

This is Greenhill's equation (1) p. 445. By it and equation (9) 
§ 6 above we have 

ajf/a^=///fifv(2)cr)(is, 

whence reducing the volume to the element d%, 

dD/dt = SV(Da) (3). 

This is Greenhiirs equation (2). Thus by equation (1) of last 
section, 

bjD^^S^a (4). 

64. To obtain Euler s equations of motion we express that 
the vector sum of the impressed forces for any volume equals the 
vector sum of the bodily forces plus the vector sum of the pressures 
on the surface. Thus 

jSSD&d% = 5SSDPd% - ffpdZ 

^ffJ(DF-Vp)d%, 

by equation (9) § 6 above. Applying this to the element ds and 
dividing by Dd% we have 

& = F-Vp/D (5), 

or by equation (1) § 62 

da/dt-SaV.a^F-VpID (6). 

This is GreenhiU's equations (4) (5) (6). 

in the text — seems to me out of harmony with the meaning attached to d in other 
branches of Mathematics. At any rate I have respectable fellow-sinners, e.g. 
Kirchhoff in his Mechaniky zweite Vorlesung, et seq.] 
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It F(pyt) = (F a. scalar) be the equation of a surface always 
containing the same particles F=0, or by equation (1) 

dF/dt-SaVF=0 (7). 

This is Greenhill's equation (7). 

65. Let us now put 

!dplD = P (8). 

This of course assumes that i) is a function of p only, which is not 
always the case, for instance in a gas where diffusion of heat is 
taking place. If P have a potential v, P = - Vv. Thus equation 
(5) becomes 

a = -V(v + P) (9), 

and equation (6) 

3(r/a^-/Sa-V.cr = -V(t; + P) (10). 

Now 

Thus equation (10) becomes 

daldt-\-2Ve<r-{-VR=:0 (11), 

where the scalar R is put for P + v — <r'72 and the vector € for 
FV(r/2. This is Greenhill's equations (8) (9) (10). 

If € = 0, i.e. Wa = 0, we have <r = V<^, whence our equation 
becomes 

V (d(f>/dt + iJ) = 0, 

so that d<t>ldt + R = H (12), 

where ZT is a function of t only. In next section we shall obtain 
in the case of an infinite fluid a generalisation of this which I 
believe has not hitherto been obtained. Here we have made the 
assumption that if Wa = at one epoch it will be so always. 
This we shall prove later. 

66. Greenhill next considers the case of steady motion. In 
this case da/dt = 0, so that equation (11) becomes 

2Vea + VE = (13), 

and therefore the surface R = const, contains both vortices and 
stream-lines and the relation dR/dn = 2qQ) sin given on p. 446 of 
GreenhilFs article is the natural interpretation of our equation. 

So far we have been going over much the same ground as 
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Hicks, but now we enter upon applications of Quaternions that I 
think have not been made before. 

67. Greenhill next considers rotating axes and finds the form 
of the equations of motion when referred to these. Let a be the 
velocity referred to them ; so that if the axes have at any time 
made the rotation q{)q~^ the real velocity will be qa-q"^. Thus*, 
as always with rotating axes, if a be any vector function of a 
particle the rate of increase of a in space is 

d + VeoUj 

where co is the angular velocity referred to the same system of 
rotating axes. Thus we see that the acceleration of a particle 
= &+ Vooa. The velocity <r = /j' + F©/}', where p is the vector 
coordinate referred to our present axes (i.e. in the notation of the 
footnote p' = qp(f% Thus our equation of motion is now 

d- + V(oa = P - p'V|)/i>, 
but now 

djdt = djdt - 8p\^ . = didt -S(a- top') p'V, 

w^hence changing />' into p we have 

dajdt -S(a-(op)V .a+ Vcoa = P - Vp/D (14), 

which is the equation Greenhill gives on p. 446. 

* This is quoted as a known result because it occurs generally in the subject of 
Rigid Dynamics. No Quaternion proof, however as far as I am aware, has been 
given. We therefore give one here. What is meant by rotating axes may be thus 
explained. — Instead of choosing as our coordinates the vectors a\ /3' ... which occur 
in any problem, we take others a, /3... such that a' = qaq~^y p' = q^q~^ .... q ( )q~^ 
may be called the integral rotation of the axes. Thus if we say that the vector 
angular velocity of these axes themselves is w we mean that the real angular 
velocity is qwf^, so that, as can be seen by putting in equation (a) below a = const., 
or as in equation (58) § 33 above, (a = 2Vq~^q. (This maybe established also by 
Tait's Quaternio7is, § 356, equation (2), from which q(aq~'^ = 2Vqq~'^ or (a = 2Vq~^q). 
Again, when we say that the rate of increase of a in space is r, we mean that 
d' = </rg~i or T = q~^ .d (qaq~^)ldt . g, or 

T=d + 2VVq-^q.a (a), 

or r = d + Fwa (6). 

This could have been proved with fewer symbols and more explanation, but the 
above seems to me the most characteristic Quaternion proof. We might have 
started with not quite so general an explanation of reference to rotating axes and 
so refrained from introducing the integral rotation, and therefore also q. 
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The Lagrangian Equations*, 

68. We now consider the history of a single particle and for 
this we require different notation. 

We consider the vector coordinate (/>) of a particle as a func- 
tion of some other vector a (say the initial value of p) and of t 

We first require the connection between ^V and pV. We shall 
drop the aflSx a and retain />, so that now not Q (V^, p^ but Q (V^, otj) 

=e(r,?). 

Now SdaV = Sdp,"^, 

but dp = — pjSdaVj, 

. • . SdaV = - SdaV^Sp^;^, 

or since da is perfectly arbitrary 

V = -V,Sp.,V (15). 

Thus operating upon equation (9) § 65 by V^/Sp^ ( ) and remember- 
ing that in that equation V must be changed to pV we get 

V^8p,&=^V(v + P) (16), 

and this is our new equation of motion (Greenhill's (1) (2) (3) 
p. 448). 

In equation (1) § 12 above let us put for p + rjy p ; and for p, ol 

Thus x« = -/Sft)V.p (17). 

Hence by equation (6d) § 14 above, 

strained vol. of el. / original vol. = S^y^s^PiPJ>jGy 
whence we see that 

DS^y,V^p^j>, = 6D, (18), 

where D^ is a constant which when a is taken as the initial value 
of p is the original density. This is the equation of continuity. 

* [Note added, 1892. At the time of writing the essay I did not notice that 
these equations are a particular case of the general equation for an elastic body 
already established (see equations (15 m) § 16 and {15n) § 16a) 



i 
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Cauchy's integrals of these eqiuitions. 

69, To obtain these integrals we require to express pV in 
terms of V. Now equation (15) of last section expresses V as a 
linear function of pV. The converse we have already seen how to 
get. In fia<;t from equation (6/0 § 3a above 

2/pV = - F/,,p^,V,V (19), 

where, with Greenhill, for brevity J is put for SVyy^p^pj)jG, 
The only function we wish to apply this to is VpV<r. We have 

= (p.Sp^'^s - PtSp/r,)S^y,V„ 

or interchanging the suflSxes 1 and 2 in the last term 

JV.Va- = p.-Srp,<r./SV,V.V. (20), 

which gives the spin at any instant in terms of our present inde- 
pendent variables a and t 

70. In order to obtain Cauchy's integral of equation (16) 
operate on it by FV ( ). Thus 

= F(V, + VJ V^Sp.a, = Sp,a, . FV.V,. 
Now 

d (Sp^a, . Wy^)ldt = Sa,<r^ . 7V,V, + Sp,&, . FV.V, . 

The first term of this last expression is zero since the sign is 
changed by interchanging the suffixes. From the last expression 
then 

^Pi^2 • ^^2^1 = constant 
= initial value of S^a,VV^^= - FVcr,, 

where a^ is the initial value of <t. Changing the suffixes and sub- 
stituting in ecjuation (20) we have 

JF,V«7 = -p,.SfVa„V., 

or giving «/ its value DJD from equation (18) nnd putting VpVa = 6, 
VVa-^ = €Qy we have 

elD = -8(eJI),)V.p (21). 

This is Oreenhiirs equations (4) (5) (6) p. 448. 

The physical interpretation of this equation is quite easy. 
Consider a small vector da drawn in the fluid initially. At the 

M. G 



82 HYDRODYNAMICS. [§ 70. 

time t this will have become dp = — SdaV . p. Thus we see that if 
a small vector ceJO^ be drawn in the fluid initially it will at the 
time t be ce/D, from which we infer that an element of a vortex 
filament will always remain an element of a vortex filament; 
or, a vortex filament or tube always remains a vortex filament 
or tube. Again we see that Te/D at any time varies directly 
with the elongation in the direction of the vortex filament so 
that Te varies as that elongation x the density, i.e. inversely as 
the cross-section of a small vortex tube at the point. This is not 
the easiest way of arriving at these results, but it is well to show 
in passing how easy of interpretation are our results. 

We see firom equation (21) that if Cq = 0, e = 0. In other words, 
if the motion have a velocity potential at one instant it will have 
one always. 



Flow, circulation^ vortex-motion, 

71. We are about to consider vortex-motion firom another 
point of view, viz. that of circulation. 

In § 12 we saw that a strain due to a small displacement 17 
could be decomposed into a pure strain followed by a rotation, the 
vector rotation being Wr)/2. If now for 77 we put the small 
vector aSt we see the propriety of calling Wa/2 the spin. This 
therefore is taken as a definition of spin. Greenhill does not take 
this (usual) course but uses the property we shall immediately 
prove concerning circulation and spin to lead to his definition. 

It is not necessary here to define flow and circulation. Putting 
a = V<^ we see that for irrotational motion the flow = —JSadp =Jd<l> 
from one point to another is the increment in the velocity potential. 
Thus for mutually reconcilable paths it is always the same. 

Taking the circulation round a closed curve in the general case, 
the curve not inclosing any singular region of the fluid, we may 
transform the line integral into a surface integral by equation (8) 
§ 6 above. Thus 

-JSadp = -JJSdXVa (22), 

so that the circulation round the curve equals twice the surface 
integral of the spin. Hence Greenhill's definition of the spin. 
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72. From equation (9) § 65 we see that 



or 



-i(/>'^'') = -t + ^ + '^'/2]^ (23). 



Therefore for a closed curve 

d (JSadp)ldt = 0, 

so that the circulation round the curve, and therefore the cor- 
responding surface integral of the spin remains always constant. 
Taking the curve round a small vortex tube we once more an'ive 
at the propositions enunciated in § 70 about vortices. 

Thus at any point of a vortex tube the strength which is 
defined as the product of the spin into the cross-section is constant 
throughout all time. Also it is the same for all points of the tube, 
for by equation (9) § 6 we have 

//SdSVo- =////SfVVdS = 

for any portion of the tube. But the only part contributing to 
the surface integral is the ends of the part of the tube considered, 
so that the strength at these two ends is the same, and is there- 
fore constant for the whole tube. 

73. These propositions are proved in the paper by Hicks 
already referred to. There is yet a third method of proof which, 
like Hicks's, is derived directly from the equation (9) of motion. 

We have from equation (1) § 62 

lv-V^ = -£f(rV.V + Vfif(rV.=V,>Sfc7,V. ...(24)*; 
therefore d VVa/dt - Wa = V ( V^fifo-.V, . o-,). 

* [Note added, 1892. When first giving this in the Mess, of Math. 1884 and 
when again putting it in the present essay, I was unaware of Prof. Tait's paper in 
Pi'oc. R. S. E. 1809—70, p. 143, where in the present notation he has for an incom- 
pressible fluid, (1) dldt = dldt - Say^ [giwen in Cartesian form], (2) o-= -^v-^pll), 
(3) Svo' = 0, (4) ^j-d{^a)ldt=-Sj^a^ .<T, (5) Vya=0, {G) d^<Tldt= - S^a^ . <t. 
Perhaps I have not interpreted Prof. Tait's notation which is but briefly described 
correctly, but (4) should apparently be Fvo--<i (vo')/^< = '^Vo'V • o"? aiid so agree 
with (5) and (6). It will be seen that the whole of §§ 512—13 of Tait's Quaternions, 
3rd ed., is contained in this essay. I cannot at present recall whether this is 

6—2 
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Now by equation (9) § 65, W& = 0, so that 

dWa/dt^rV^tT,S^^a, (25), 

.-. d{p + xWa)ldt = <T + xVVa- + xW^a^^^i, 

Now a-,SV,<T,V^= Va^V^SV^a^ + FV, 7^0-^(7,+ W^a^^a^, 
but T^i^jSo-jO-, = 0, for an interchange of suffixes leads to a change 
of sign. Therefore 

FVjO-^VjO-j = FVo-fifVo- - SVaV . <r. 
Put now X = c/D where c is some small constant. Thus 

. xlc = -i>II>^-'8V<rlD, 
and we get d(p + cVVa.lD)dt = o- - cS^aV .a./D (26). 

This shews that p =p-\' c WajD is the variable vector of a 
particle, for the equation asserts that the time-flux of p' = the 
velocity at p\ Thus again we get the laws of vortices given in 
§70. 

Irrotational Motion, 

74. We use this heading merely to connect what follows with 
what Greenhill has under the same on p. 449. It is not very 
appropriate here. 

For irrotational motion we have seen that we may put 

o- = V<^ (27). 

If the fluid be incompressible we have further 

= SV(r=V^<^ (28). 

Let T be the kinetic energy of this liquid. Thus 
- 27= DJJf{V<l>yd% = DJJ<l>8d%V<l> - DfJJ(l>V^d% 
by equation (9) § 6 above. Thus by equation (28) 

T='-iDff<t>SdlV<f> (29). 

Hence we see that if we have a vector a which satisfies the 
equations FVo- = and 8Va = throughout any singly-connected 
space, or more generally satisfies the second of these equations 

owing wholly or in part to my being indebted to other old papers of Prof. Tait, or 
whether in writing his 3rd ed. he has arrived independently at the same treatment, 
but am inclined to the latter belief.] 
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throughout any space, and also has zero circulation round any 
closed curve in that space*; and further satisfies the equation 
SdS,a- = at the boundary ; then the function 

JJJ<T'd% = 0, 

or <r = at all points of the space, for each element of the integral 
being essentially negative must be zero. 

75. The following important theorem now follows : If there 
are given ; at every point of any region the convergence SV<r and 
the spin VVa, at every point of the boundary the normal velocity 
(and therefore SdX<r), and the circulation round every cavity which 
increases the cyclomatic number of the region ; then the motion 
(as given by a- at every point) if possible, is unique. For let a be 
one possible velocity, and if possible <r+T another. Thus t 
satisfies all the conditions that <r does at the end of last section, 
and therefore t = at every point or <r is unique. 



Motion of a solid through a liquid. 

76. We assume that there is no circulation of the liquid for 
any cycje ; in other words, that if all the solids be brought to rest, 
so will be also the liquid. 

We shall take axes of reference fixed in the (one) moving solid. 
In the foot-note to § 67 is explained how the rotation of these axes 
is taken account of. The effect of the translation of the origin 
will be easily found. 

Let <r, (0 be the linear and angular velocities respectively of 
the moving solid. Let us put 

<l> =^ - Sayjr - S(OX (30), 

where <f> is the required velocity potential of the liquid, and -^ and 
^ are two vector functions of the position of a point independent 
of a and w. Let us see whether yjr and x oan be found so as to 
satisfy these conditions, a and co are of course assumed as quite 
arbitrary. 

The conditions are first the equation (28) of continuity 

V2<^ = (28), 

* To insure that <f> is single valued and so, that the part of the surface integral 
of equation (29) due to " barriers " is zero. 
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which gives 

V2^ = 0, V^x = (31), 

and second, the equality of the normal velocity —8UdZV(f} of the 
liquid at any point of the boundary with that of the boundary at 
the same point. Thus for a fixed boundary 

So- (SdXV .ylr) + Sco (SdXV . x) = 0, 

whence on account of the arbitrariness of <r and © 

£fd2V.i|r = 5d2V.x=0 (32). 

For a moving boundary again we have 

SdXV^ = Sd% (<r + Vcop), 

or So- (dS + Sdl.V . f) + StoiVpdX + Sd^V . %) = 0, 

which gives 

dX+SdV7 .yfr = (33), 

VpdX + SdXV.x^"^ (34). 

Now it is well known that yjr and ;)^ can be determined to satisfy 
all these conditions. In fact oo being a coordinate of either -^ or ^ 
these conditions amount to: — V'^x = throughout the space, and 
SdltVx = given value at the boundary. 

77. We do not propose to find i|r and x ^^ ^^J particular 
case. We leave p. 454 of Greenhill's article and go on to p. 455, 
i.e. we proceed to find the equations of motion of the solid For 
this purpose we require the kinetic energy of the system. Calling 
this T we shall have 

T= - S(r'Zal2-' Sco^(r - ScoClcol^ (35), 

where 2, ^ and 4> are linear vector functions and the first two are 
self-conjugate. This is the most general* quadratic function of 
ft) and 0-. It involves 21 independent constants, six in 2, six in Xi 
and nine in <I>. When yjr and 'x^ are known 2, ^ and H are all 
known. We will obtain the expressions for them, and for simplicity 
will assume that the origin is at the centre of gravity of the 
moving solid. Let M be the mass of this last and fio) (where, as is 
well-known, /^ is a self-conjugate linear vector function) the 
moment of momentum. Thus by equation (29) 

2T= - DJJ^SdXVcp - Ma' - ScjfjLco. 

* For let A (o-, o") + i^ (w, <t) + G (w, w) be this general function Aj i?, and C being 
scalar functions each linear in each of its constituents, and let 
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Putting 8dX^<f> = 8d% {a + Vtop) at the moving boundary and 
zero at the fixed, 

where the surface integral must be taken only over the moving 
boundary. Thus noting that S and 12 are self-conjugate 

X\ = M\-^DJJ(ylrS\d% + dl8\ylr) (36), 

ax = /^X - ^DJJixSXpdl. 4- VpdlSXx) (37), 

- ;Sf\4)V = ^DJJiSXxS^'dX + SXpdXSX'f) (38). 

These surface integrals can be simplified, for in each of these 
equations the first surface integral equals the second. . In equation 
(36) in order to prove this we have merely to put for d%, — 8dl,V .yjr 
(equation (33)) when we shall find by equation (9) § 6 above 
(since we may now suppose the integrations to extend over the 
whole boundary), that 

Similarly for equation (37). Again in equation (38), 

//iSfXx^'dS = -JJJSVy,S\xiS)^'i^2d^ = JJSXpdXSX'ylr. 
Thus finally for 2, 12 and 4) 

XX = MX- DJJfSXdX = MX- DJJdXSXylr (39), 

nx = /jX-' DJfxSXpdX =fiX- DfJVpdXSXx . . .(40), 
_ Sx^X' = DJfSXxSX'dt = DJJSXpdXSx'yIr (41). 

This last may be put in the following four forms 

^'X = -DJIdtSXx = -DJfypSXpd% (42), 

<E>X = - DJfxSXdl. = - DJJVpdXSXyJr (43). 

Thus assuming that yjr and x ^^^ determined we have found T 
as a quadratic function of a and cj. 

78. If P, G be the linear and angular impulses of the system 
respectively our equations of motion are 

P+ Fft)P = F, 
G+FcoG-M, 

by the footnote to § 67 above. Here F and M are the external 
force and couple applied to the body. Now at the instant under 
consideration P = ^VT, G-=^VT. But if the origin had been at 
the point — p, P would still be ^VT whereas G would be ^^T+ VpJ?, 
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Thus diflferentiation with regard to t does not aflfect the form of P 
but does that of G. In fact using the last stated values of P and 
G along with the last two equations and eventually putting 
/) = 0-, /o = we get 

d^VT/dt+ Fo)^Vr=P (44), 

dST/dt+V(oJ7T-{- VaJ^T^JH (45). 

Now from equation (35) we see that 

,Vr=2<r + a>'o, (46), 

«Vr=^(r + Ila) (47), 

whence from equations (44) and (45) 

X<r + ^'o) + FwSo- + Fa)4)'(» = P (48), 

^(T + Ho) + Fo-So- + ( F<r* © + Fa)*(r) + Fwflo) = M. . .(49). 

Thus we see that with Quaternion notation even the general 
equations of motion are not too complicated to write down con- 
veniently. 

We now leave Greenhiirs article and proceed to certain 
theorems not contained therein. The Cartesian treatment of 
these subjects will be found iu Lamb's treatise on The Motion of 
Fluids, Chapters Vi. and IX., which are headed Vortex Motion and 
Viscosity respectively. 

The velocity in terms of the convergences and spins, 

79. We have seen in § 75 that if the spin and convergence 
are given for each point of a bounded fluid and the normal 
velocity at each point of the boundary, there is if any, but one 
possible motion. We shall see directly that this unique value 
always exists. 

At present we observe that we cannot lind a motion giving an 
assigned spin and convergence for each point and any assigned 
velocity for each point of the boundary. If however with such 
data a motion be possible we can find the velocity at any point 
explicitly. By equation (19) § 10 above we have 

47ra = - VJJJVuad%, 

but by equation (9) § 6 

Jff^ii(rd% =JJudX(T -.JfJiiVad% 

.-. ^7r(r^JJ[V]VvdXa-JJJ[V]VuVad% (50), 



§ 80.] HYDRODYNAMICS. 89 

where the square brackets indicate that we may or may not 
retain the Fat our convenience. This equation may be put 

iira- =// [F] (VudXa - udXVa) -^fffuV^adS (51). 

Both of these equations solve the question now proposed. 

If the fluid be considered infinite the surface integral of 
equation (50) vanishes if a- converges to zero at infinity and also 
that of equation (51) if in addition the spin and convergence 
converge to a quantity infinitely small compared with the re- 
ciprocal of the distance of the surface. 

80. We may now consider the case when V<r (spin and con- 
vergence) is given at all points and Sd^a- (normal velocity) at the 
boundary. It must be observed that the given value of the spin 
must be distributed in a solenoidal manner for 

;SfVFV<r = 0. 

Whenever the quaternion 

^iTq^-SJSVuVad%^VjjjuVad% (52), 

is a vector, we see that all the conditions except the boundary 
ones are satisfied by putting <r = g'. Let us then see when q 
reduces to a vector. By equation (9) § 6 above 

^TrSq = - JJuSdtVa + JfJuS^'ad%, 

Thus 5 is a vector if this surface integral vanish. The surface 
integral vanishes if all the vortices form closed curves within the 
space. If they do not we must extend the space and make them 
form closed curves outside the original space. Extending the 
volume integral accordingly, we may put 

47r(r'==-JJfVuVad% (53), 

and now Va-' = Vcr. Suppose then that 

a = a'+<r'' (54). 

We thus get V<r'^ = and may therefore put 

(7" = V0 (56), 

where <^ satisfies the equations 

V^(t> = (57), 

and 8d^V(l> = SdX (<r — a) = known quantity (58). 

Now it is well known that <^ can be determined so as to satisfy 
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these two equations. Therefore the problem under discussion 
always admits of solution. 

The above is equivalent to Lamb's §§ 128 — 131. His § 130 is 
the natural interpretation of the equation 

4 TTo- = - /// WuVadS. 

His § 132 is seen at once from equation (50) above. For in the 
case he considers we, in accordance with § 7 above, take each side 
of the surface of discontinuity as a part of the boundary. Now 
[SdXa]a^i, = so that for this part of the boundary we can leave 
out the part Sdlto- and we get 

^Tra = JJVVurdta-JJJVVuV(rd% (59), 

so that if we regard [FdScr]^^^ as — 2 x an element of a vortex, 
we get the same law for these vortex sheets as for the vortices in 
the rest of the fluid. 

81. The velocity potential due to a single vortex filament of 
strength dO is at once obtained by putting in equation (50) for 
FVcrdS, 2d0dp. Thus calling cr' the part of the velocity due to 
the filament 

27ro-' = - dejwudp = - dejju.v^^vdtv^ 

by equation (8) § 6 above. Now 

V^.VdtV^ = - dtV,' + V.fifdSV,. 

Hence, since V^it = for all points not on the filament 

27r(7' = V (dejfSdXVu) (60). 

Thus the velocity potential = (27r)"^ x the strength x the solid 
angle subtended by the filament at the point considered. 

Kinetic Energy. 

82. Let us put 

4<7rq=fJJaVad% (61), 

so that (T = Vq (62). 

We may now find expressions for T, the kinetic energy, in one 
or two interesting forms. We have 

2T=-JffDa'd%, 

whence 2T = - JJJDS(7Vqd% = -fJDS<rdlq +JJfD^8ay,qd%. 
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We assume that the fluid extends to infinity and that the 
vortices and convergences are all at a finite distance, so that at 
infinity a is of order \\R^ and 3 of order l/E. Thus the surface 
integral vanishes and we have 

2T^SSJD,S<T^V^qd% (63). 

or, putting in the value of q from equation (61), 

T = (Stt)-' //////MA'S<r.V.V,«r,dMS, • -(64). 

These are Lamb's equations (28) and (29), p. 160, generalised. 

Similarly his equation (30) generalised is 

2r = fffD^8pay,<y,d% (65), 

for JffD,Spa,V,<T,d% = SJDSp<7dta- - fffD8^<r^ad%, 

by equation (9) § 6 above. But 

?<T?=25S<r?-r=<7=<7, 
and the surface inteojral vanishes as before. 



Viscosity. 

83. To consider viscosity the assumption is made that the 
shearing stress which causes it is /^ x the rate of shear of the 
moving fluid, fi is assumed independent of the velocity and 
experiment seems to shew that it is also independent of the 
density. This last we assume though the variation with density 
can be easily treated. 

Consider a general strain %a). Since 

^Xf^X^ = ^x"' ^Z^^- 
Putting fi and v for any two vectors perpendicular to co we see 
that the normal U(o to any interface becomes Ux~^(d by the strain. 
Hence the interface &> experiences a shear (strain) which equals 
the resolved part of x^^ perpendicular to the vector x~^(o equals 
resolved part of (% — x~^) Uo) perpendicular to x'^co. Now when 
X is the strain function due to a small displacement adt, by § 12 

above 

XO) = (o — ScoV . (rdt (66), 

* Kelland and Tail's Quaternions^ chap. x. equation {n). We have already used 
a particular case of this in § 15, above. 
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whence X '*^ = (o-^V^Scja^dt (67), 

and we see that the shear is the resolved part of 

perpendicular to w, i.e. parallel to the interface o). From this we 
see by our assumption concerning viscosity that any element of 
the fluid is subject to a stress <^ given by 

jB being a scalar. Now we define the pressure p by putting 

.-. (l>(o = -pQ} + ^fjLwS^a- - fjL (S(oV . a + V^Scoa^), . .(68). 

The equation of motion is 

i)(r = DP + <^A, 
or 

D (da/dt - SaV . 0-) = Do- = DP - V^ - fiVS^a/S - /^vV. . .(69). 

If fjL be not as stated independent of D this equation must be 
made to contain certain space fluxes of /jl which are quite easy to 
insert. 

84. In § 179 of Lamb's treatise he considers the dissipati(m 
function due to the viscosity of a fluid in a way that seems to me 
misleading if not wrong. It appears as if he should add to the 
first expression of that section 

itdpxx/dx + vdpxy/dx + wdpxzjdx. 

He refers to Stokes, Camh. Trans, vol. ix. p. 58. Let us give the 
quaternion treatment of Stokes's method. 

If T is the kinetic energy of any portion of the fluid 

Thus •.• d{Dd%)ldt^Q, 

f = -JJJD8a&d%, 

but Do- = DF + <^iV,, 

so that r = - JJJDSaFdS - JJJS(r(l>y ^d% 

or f = - JJJDSaFd% - JJS(r(t>d^ + IIIS(T,<t>V^d% (70). 

If now (f) be expressed in terms of p and /jl we have t depend- 
ing on F, p and /a. The part of T depending on F and p repre- 
sents energy stored up as potential energy of position and potential 
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energy of strain respectively, but the part depending on fi repre- 
sents a loss of energy to the system we are considering the energy 
being converted into heat. 

Thus putting <^ =^ + w we have by equation (68) 

= f fjLcoS^f^ + 2fifco J ^ ^' 

where yjr is given by equation (75) below. 

Thus we see that the rate of loss of energy is 

//S«rW2 - fffSa,wV^d% (72). 

The surface integral is the work done by viscosity against the 
moving fluid at the boundary and the volume integral is con- 
sidered due to the work done against the straining of the fluid. 
Thus we put 

F^^Sct^t^tV^ (73), 

and call F the " dissipation function." 

By equation (6) § 3 above 

i^=-£f'^?«r? (74), 

where ylr is the rate of pure strain of the fluid, i.e. 

2yft(o = ^S(oV,(r-V^S(0(r^ (75). 

Again by equation (18) § 18 above because F is quadratic in yjr 

F=-s>^Uani^ (76), 

so that from equation (74) we have, by § 4 above, 

^^^am (77). 

Substituting for ^ from equation (71) in equation (74) 

F 1 ^ {SKHf - 1,,^m- (78), 

which gives F in terms of y^. 



SECTION VI. 
The Vortex- Atom Theory, 

85. If Quaternions can give valuable hints or indicate a 
promising method of dealing with the highly interesting mathe- 
matical theory of Vortex-Atoms, I think this alone ought to be 
suflBcient defence of its claims to be within the range of practical 
methods of investigation. 

In what follows I think I may be said to have indicated a 
hopeful path to follow in order to test to some extent the sound- 
ness of this theory. 

Statement of Sir Wm, Thomsons and Prof, Hicks^s theories, 

86. Sir Wm. Thomson's theory is so well-known that it is not 
necessary to state it in detail. Matter is some differentiation of 
space which can vary its position carrying with it so to speak 
certain phenomena, some of which admit of definite quantitative 
measurement. Perhaps the most important of these phenomena 
is what is called mass. Now, says in effect Sir Wm. Thomson, if 
we suppose all space filled with an incompressible perfect fluid the 
vortices in it are just such differentiations of space. They also 
carry about with them definite quantitative phenomena. Can we 
prove that these hypothetical vortices would act upon one another 
as do the atoms of matter, the laws of whose action are contained 
in the various Sciences, e.g. Physics, Chemistry and Physiology ? 
The problem in its first stages at any rate is a mathematical one, 
but during the many years it has been before the mathematical 
world very little progress has been made with it. 

Hicks's extension of this theory is perhaps not so well known, 
but it seems to me quite as interesting and more likely to tally 
with the known phenomena of matter. He enunciates his theory 
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in the Proceedings of the Cambridge Philosophical Society y vol. ill. 
p. 276. It differs from Thomson's simply in assuming that the 
fluid does not quite fill space — that there are in it bubbles*. 
These bubbles will find their way to. where the pressure is least, 
i.e. speaking generally to the centre of some at least of the vor- 
tices. Thus we have another source of differentiation of space 
and general considerations seem to point to these differentiations 
being the true atoms, though of course "atom" is no longer a 
descriptive term. 



General considerations concerning these theories, 

87. In Max weirs article Atom in the Encyc, Brit, p. 45, he 
says — "One of the first if not the very first desideratum in a 
complete theory of matter is to explain first mass and second 

gravitation In Thomson's theory the mass of bodies requires 

explanation. We have to explain the inertia of what is only a 
mode of motion and inertia is a property of matter, not of modes 
of motion. It is true that a vortex ring at any given instant has 
a definite momentum and a definite energy, but to shew that 
bodies built up of vortex rings would have such momentum and 
energy as we know them to have is in the present state of the 
theory a very difficult task. 

" It may seem hard to say of an infant theory that it is bound 
to explain gravitation." 

Now as Hicks tells us what induced him to give his theory 
was the promise it gave of explaining gravitation. But I believe 
nowhere does he point out the still more important result that 
probably on his theory we can explain inertia. 

The statement of the principal property of inertia put scien- 
tifically is that the motion of the centre of gravity of any two 
bodies approximates more and more nearly to uniform velocity in 
a straight line, the more nearly they are isolated from external 
influence. But this property is probably true of Hicks's bubbles. 
The centre of gravity of any portion of the fluid containing certain 
bubbles (1), (2) ... {n) will if approximately isolated from all the 

* By ** bubbles " of course I do not mean spaces occupied by another kind of 
fluid of smaller density but actual vacua. Thus a bubble may start into existence 
where none previously existed, or again a bubble may completely disappear. 
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rest of the fluid move approximately in a straight line, but this 
amounts to saying that the centre of volume of the n bubbles will 
also move uniformly in a straight line if the centre of the whole 
volume (bubbles and fluid) considered, move similarly. These 
conditions are not evidently true, but I think they are probably 
so when we consider groups of large numbers of bubbles. 

Whether this theory explains gravitation is one of the principal 
questions to be considered in the following first trial. 



Description of the method here adopted, 

88. Finding it practically impossible to consider the real 
problem of a number of bubbles in the liquid, I consider the fluid 
continuous and not incompressible. Let us assume that 

D = ianh(p/c) (1), 

where c is some very small constant which in the limit = 0. For 

ordinary values of p, D very nearly = 1. It is only when p 

becomes comparable with c that D varies. When p = 0, -D = 0. 

Also 

P=/(dp/i)) = clogsinh(p/c) (2). 

For ordinary values oip then, P=p but when p becomes compar- 
able with c, P varies and when jp = 0, P = — oo . If we assume 
then for the greater part of the fluid that p is large compared 
with the small quantity c we see that the fluid will have almost 
exactly the same properties as a liquid containing bubbles. We 
may now apply the equations of motion of §§ 64, 65. 

Now we know the velocity at any point of an infinite fluid in 
terms of the spins and convergences at all points. From this we 
may deduce the acceleration in terms of the spins, convergences 
and time-fluxes. It so happens that by the equations in § 73 we 
can get rid of the time-fluxes of the spins. This greatly simplifies 
the further discussion of the problem. The equation that we thus 
deduce forms the starting point of our investigation, for the phe- 
nomena of gravitation, electro-magnetism, stress, &c., are exhibited 
and measured by the acceleration in bodies due to their relative 
positions. 

The equation we obtain at once gives a generalisation of the 
integral, equation (12) § 65 above. 

At the end of this section I give an equation which is rather 
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complicated but which promises to enable us to deal more 
rigorously with our problem than I profess to have done below. 

We proceed to the investigation just indicated. 

Acceleration in terms of the convergences, their time- fluxes y and 

the spins, 

89. Let us put for the convergence and twice the spin m and 
T respectively, so that 

SV(r = m (3), 

FVcr = T (4). 

We have seen that in equation (.50) § 79 we may neglect the 
surface integral. Thus 

4^7ra- = VJJJu(m + T)d% (5), 

.-. 4^7rda/dt = VJJJud (m + T)/dt . d%. 

By §73 T=TSV(T-8TV.a, 

but by equation (1) § 62 

dr/dt = f + SaV . T, 
.-. dr/dt = SaV . T + rSVa- - SrV . a, 

or because SVt = 

dr/dt = TSAa-(rS^T=WV(TT (6) 

.-. 4^7rd(rldt = S/JfJu{VS;V(jr-\-dmldt)d% (7). 

This equation is not in a convenient form for our purpose, for 
T and m may be and most probably will be discontinuous, so that 
it is advisable to get rid of their derivatives. Moreover it is 
advisable to allow the time-flux of m to appear only under the 
form d (md%)/dt because in case of pulsations the time integral of 
this expression will be zero, whereas we can predicate no such 
thing of (dm/dt)d%, or indeed of 7hds. Let us first then consider 
the first term in equation (7), viz. VJJJuV^Va-rdS. Using equation 
(9) § 6, and neglecting the surface integral at infinity, this 
becomes 

- VJJJVVu VardS = V WJffu Vard% [§ 9] 

= ^7//^' y(rrd% - VSS/JJJu VardS 

= 47r Var + VJJJSarVudS, 

by equation (19) § 10 above. 

M. 7 
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Remembering now (§65 above), that 

&=daldt-VaT-V{a^)l2, 
we see by equation (7) that 

47r(7 = - 27rV . 0-2 + VflfSarVudS + Vfjftt (dm/dt) d% . . .(8), 
whence by equation (9) § 65 above 

P + v - <ry2 + (4iTry'JJJ(SaTVu + vdm/dt)d% = H .. .(9), 

where iT is a function of the time only. This is a generalisation 
of equation (12) § 65, for assuming that t = we know by § 79 
that d(l>ldt = {4i'rry'fjfu (dm/dt) d%. 

90. This integral equation may be put into several different 
forms by means of equation (9) § 6 above. The form that is 
useful to us is the one in which instead of dm/dt we have d{md%)/dt 
as we have already seen. Now d (d%)/dt = — md%. Therefore 

d {md%)/dt = (m — m*) d% 

= (dm/dt - Sa-Vm - m*) d%. 

Substituting from this for dm/dt, and then transforming by 
equation (9) § 6 so as to get rid of the space variations of m in- 
volved in SaVm we get 

JJJu (dm/dt) d% = JJJvd (md%)/dt + /// (um^ - mS<rVu - um"") d% 

= Sjjud (md%)/dt - JJJmSaVud%. 

Thus from equations (8) and (9) 

47ra=-27rV . a'-\-S/JJf8S/u(VaT-ma)d%+S/JfJud(md%)/dt.,,('^0X 

and 

P-{-v-a'/2-\-(4^7r)-'JJJ[d%SVu(V(rT'-ma-)-{-ud(mdS)/dt}=^H...(ll), 

Sir Wm. Thomson's Theory, 

91. We are now in a position to examine the two theories. 
We first take Thomson's, which is considerably the simpler, and 
which therefore serves as an introduction to the other. We have 
then m = 0. Thus equations (10) and (11) become 

47r6- = - 27rV . crV V/Z/fifo-TVitcis (12), 

p/D - (772 + {^7r)-'ffJSaTVud% = H (13), 

for in the present theories we may put ?; = 0. 

We shall consider the two terms on the right of equation (12) 
separately. The second term gives then an apparent force per 
unit mass due to a potential 

-(^irT'jjjSarV^idS (14). 
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Comparing this with equation (33) § 46 above we see that this 
potential is the same as that of a magnetic system given by 

47rl= Var (15). 

Now this magnetism is zero where t is zero. In other words 
it is only present where the vortex-atoms are, and therefore it 
cannot be so distributed as to give an apparent force of gravitation, 
for taking the view of magnetic matter expressed in equation (34) 
§ 46 we see that there must, in every complete vortex-atom, be 
a sum of magnetic matter exactly = 0. Nor again is it likely to 
explain the phenomena of permanent magnets, because, assuming 
that for any given small space including many vortex-atoms 
JJJV<rTd% is not zero, the apparent force produced will affect all 
other parts of space independently of whether this same integral 
for them is not or is zero. But to explain the phenomena of 
permanent magnets we must assume that the effect takes place 
only on portions of space where there is positive magnetic matter. 
This term then gives us no phenomena analogous to physical 
phenomena. As a matter of fact it probably has no visible effect 
on large groups of vortices, for there is no reason to suppose that 
the vector Vcrr is distributed otherwise than at random. 

92. Let us now consider the other term in equation (12), and 
neglecting the term already considered, put 

o- = -V.<r72 (16). 

The phenomena resulting from this are the same as would 
follow from a stress in a medium, the stress being an equal tension 
in all directions = — a^/2. Now comparing equation (62) § 58 with 
equation (4) § 89, we see that a depends on T/47r in exactly the 
same way as does H on C. The question then arises — is the stress 
we are now considering equivalent in its mechanical effects upon 
the vortex-atoms to the stress given by equation (68) § 60 above, 
which explains the mechanical effect of one current on another ? 
We saw in § 60 that this stress is a tension — J-V/Sir along the 
vector H and an equal pressure in all directions at right angles. 
The effects then would be the same only if a be at right angles to 
the surface of our atom. But this is obviously not in general the 
case. From this analogy we can see however what approximately 
will happen to our atom. For instead of a being at right angles 
to the surface it is in all probability very nearly tangential. 
Assuming that it is actually tangential we see that at the surface 
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of the atom we have a tension exactly corresponding to the 
pressure which in the electric analogue will be exerted on this 
surface. In other words, the atoms will act on each other very 
approximately in what may be called a converse way to the small 
circuits in the electric analogue; i.e. where, in the electric analogue 
there is an attraction, in the hydrodynamic case there will be an 
apparent repulsion, and vice vers^L. 

Now each vortex-atom forms a small circuit and therefore acts 
in the converse way to a small magnet. In other words, each atom 
acts upon each other atoni as if it were charged with attracting 
magnetic matter. Thus we see that if we could suppose certain 
extra atomic vortices to exist and to be disposed throughout space 
in what at present must be considered quite an artificial manner 
with reference to the atomic vortices we could rear up a fabric 
which would explain gravitation. This conception however is of 
very little use for our present purpose. 

From these considerations I think we have every reason to 
believe that Thomson's theory in its native simplicity does not 
promise to lead us to the physical phenomena of matter. We 
pass on therefore to Hicks's. 



Prof. Hicks's Theory, 

93. Hicks in his theory, as I understand him, assumes that 
the bubbles always remain associated with the same particles of 
the fluid. This of course is probably not the case. By reason of 
the variation of the pressure with the time it is probable that 
evanescent bubbles start into existence and disappear at various 
parts of the fluid. This requires some few preliminary remarks. 

The particles of the fluid with which bubbles are permanently 
(i.e. throughout the greater part of each small but not infinitely 
small interval of time) associated are those where the intensity of 
spin is greatest. If the intensity of spin is quite various at 
different points we shall thus have vortices where there is 
generally no bubble, extra material vortices in fact. We must 
suppose these distributed quite at random till the more exact 
mathematical treatment of our problem leads us to suppose other- 
wise. Now evanescent bubbles will occur rather in these vortices 
than in parts of the fluid where is no vortex at all (if we may 
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suppose such parts to exist), and of course they will occur more 
readily in stronger than in weaker vortices. At the present stage 
of the theory then we may suppose evanescent bubbles to occur in 
all parts of the fluid. As a first approximation to the considera- 
tion of the effect of these bubbles we may assume a part rn! of m 
to be continuously distributed through space. Putting 

m = m! + M (17), 

we must suppose M to be present only at the material bubbles 
where it is probably discontinuous, whereas rnf is continuous 
throughout both the material vortices and the rest of the fluid. 

Now when on account of variation of pressure m' and M are 
affected — is it probable that in the neighbourhood of a permanent 
bubble m' and M are of the same or opposite signs ? To answer 
this question observe that what we call m continuously distributed 
is really a series of discontinuous values of m scattered at random 
through space so that m' is probably very small. A decrement of 
pressure will cause an increment of evanescent bubbles, i.e. a 
decrement in in. An increment of the permanent bubbles will 
also take place the magnitude of which by what has been said 
concerning rnl will not be by any means accounted for by the 
decrement in m. There will therefore also be a decrement in M, 
Similarly for an increment in the pressure. If and m' may there- 
fore be assumed to be of the same sign. 

After noticing that m is continuous and therefore that there is 
no objection to introducing its space variations we are furnished 
with all the materials necessary for discussing our problem. The 
equation we shall use is (10) of § 90 above. We divide its discus- 
sion into two parts as follows. 



Consideration of all the terms except — V(a^)/2, 

94. The reasons that we have already seen in § 91 for 
neglecting VJJJSo-TVud% still hold good so we put this aside. 
This is not the case with —VfJJmSaVud% but we can neglect 
^JJJud {md%)/dt for the average value of d (md%)ldt for any particle 
is zero. The only term to consider then is —VJJJmSaVudS. 
Putting as in equation (17) m = m' + M we have 

- VJJJMSaVitdS - Vf(fniSaVud%, 
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The first term of this can be neglected for the same reasons as 
for neglecting that containing F<7t. Applying equation (9) § 6 to 
the last term and neglecting the surface integral as usual we get 

VJJJuS^ (mV) d% = VJJJu {m' (m -hM)-h /SVVm'jdg. . .(18). 

The last term can probably be neglected though we cannot give 
such good reasons as for the other terms we have neglected. At 
any rate in places not near permanent bubbles SaVm is as likely 
to be positive as negative and vice versA, so that such portions of 
space will on the whole produce no effect on the permanent 
bubbles. If SaVmf contributes anything for parts of space in the 
neighbourhood of permanent bubbles we must be content at 
present with the assumption either that the contributiou is in 
general positive or that if it be negative it is not sufl&cient to 
cancel the effect of the positive term m'M. Remembering that 
m is small compared with M we are left with the positive term 
mfAI. This as can be easily seen from the form of equation (18) 
leads to an apparent law of gravitation for our perjnanent 
bubbles. 

95. The gravitational mass which we must on this supposition 
assign to each permanent bubble varies as the average value of 
mM for that bubble*. Now we saw in § 87 that the probable 
measure of mass of a permanent bubble was proportional to its 
average size. Do these two results agree ? I cannot say, but even 
if they do not these considerations would still explain the motions 
of the solar system, but if the sun and Jupiter (say) were to 
collide their subsequent motion would not be that due to the 
collision of two bodies the ratio of whose masses is that which is 
accepted as the ratio of the sun's and Jupiter's. As a matter 
of fact however I should imagine that the average value of mM 
for a permanent bubble is proportional to its average volume and 
this simply as a consequence of the reasoning in § 87 above. 

A conclusion at any rate to be drawn from the above is that 
there is a presumption in favour of Hicks's theory explaining 
gravitation. 

* There is one important difference to be noticed between this and Hicks's 
explanation of gravitation. His depends on the synchronous pulsations of distant 
vortices. I do not wish to imply that I do not believe in the existence of such 
synchronous pulsations, but by the above we see that gravitation can probably 
be explained independently of them. As a matter of fact such synchronous 
pulsations probably actually occur on account of the variation of H with the time. 
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Consideration of the term — V . <7'/2. 

96. In considering this term we adopt the method of § 92 
and consider an electric analogue. The analogue is an electro- 
magnetic field for which in the notation of §§ 46 to 60 above at 
every point. 

H = o- (20). 

For this field we have at once 

47rC= FVH = T (21). 

The distribution of the magnetism in the field is somewhat arbi- 
trary, but in the notation of equations (61) and (62) § 82 it will be 
found that everything is satisfied by putting 

4^7rl = -VSq (22). 

This gives as it should S^ (H + 47rl) = 0, which is in fact the only 
equation it is necessary to satisfy. We have further 

B = H + 47rl = o- - VS^r = V F^. 
so that A= Vq (23). 

Thus all the important vectors in the analogue are determined. 
It remains to compare the mechanical effects of the analogue with 
the term - V . 0-72. 

I, it must be observed, is not confined to the bubbles, but is 
distributed throughout space. 

97. If we now assume that bubbles have not always existed 
in the positions which we call permanent, there cannot at the 
surface of the bubbles be any circulation round them. This makes 
the velocity at the surface almost normal to it, so that the stress 
given in equation (6G) § 60 reduces to a tension — SH (2B — H)/87r 
over the surface, i.e. we have a pressure 

Now on account of the absence of circulation Wq is very small 
and may therefore be neglected. Thus we get a pressure — (t^/Stt, 
and are thus led once more to a "converse'* of the analogue. 
This at once* leads to another reason for the law of gravitation if 

* [Note added, 1892. Because the density of attracting magnetic matter of the 
analogue = - *SvH/47r = - 7«/47r.] 
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the pulsations are synchronous. This we have already seen to be 
probable. 

The present consideration of the subject is merely to point to 
a method of investigating the theory of vortex-atoms. I therefore 
leave the subject here, not attempting to force the phenomena we 
have been considering to tally with the known phenomena of 
electricity and magnetism. Nevertheless I may say that the 
prospect of discussing these things by means of the present subject 
can scarcely be considered as distant after what has gone before. 

To sum up, this first application of the method leads to a pre- 
sumption in favour of Hicks's theory leading to an explanation of 
both the important properties of matter — inertia and the law of 
gravitation — and there is also reason from it to hope that the 
phenomena of electro-magnetism are not unlikely to receive an 
explanation. Thomson's theory on the other hand would seem to 
fail in the first two at any rate of those endeavours. 

98. We close the essay with the fulfilment of the promise 
made towards the end of § 88. In that section it will be remem- 
bered we considered a hypothetical fluid for which D = tanh (p/c), 
and made this do duty for a liquid containing bubbles. Strictly 
speaking our liquid is bounded at the bubbles and therefore as a 
bounded liquid should it be treated. For such a liquid we require 
an equation corresponding to equation (10) § 90, and if possible 
equation (11) also. This last I have been unable to obtain, and I 
am not sure that to solve the problem explicitly is possible. 

Our problem only deals with an incompressible fluid, but as 
the removal of this restriction does not greatly complicate the 
work we will consider the general case of a bounded compressible 
fluid. We have 

47ro- = V\[[fuad% = - V JJJVuad% 

= V ///^^ (t + m) dS - V JJud^a, 
whence 

4^7rd(rldt = V JJJd {u(t-\- m) d%}ldt - V Jjd {udta-)/dt 

The justification of using d/dt on the left and d/dt under the inte- 
gral sign will appear if the increment (dcr/dt) dt in o- at a given 
point in the time dt be considered. It will be observed that the 
meaning here to be attached to u will be —SaVu, as in the 
differentiation d/dt with regard to the time the origin of tt is 
assumed to be fixed. 
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I shall now merely indicate the method of procedure. By the 
method exhibited in § 89 we can prove that 

d (uTd%)/dt = - WuVa'Td% - uoStMs 
and that 

d (umd%)/dt = — m8(rVud% + vd (mds)/dt. 

From this we can deduce that 

4m'da/dt = VrVJJJuV<rrd% + V /// {- mSoVudS + ud (md%)ldt] 
- VjjuaSrdX - Vjjud {dto-)/dt + V fJdilaSaVu, 

and from this again we get 

&=zVw + Vw' + Vv (24), 

where w and w' are scalars and v a vector given by 

4iTrw = - 27ro-' +/// [d%SVu ( Var - ma) + ud (md%)/dt} . . .(25), 

so that w is in fact the H — v — Poi equation (11) § 90 

^7rw'=JJ(8dta8(rVu-ud8dta/dt) (26), 

47ri; =//(-^ uaSrdt + VdlaSaVu - udVdXajdt) . . .(27). 

This last equation may be put in what for our purposes is the 
more convenient form 

47ri; = - // {uaSTdt -\-d{u 'Vdl(T)ldt] (28). 

Again it may be put in a form free from djdt', for FdSo- = 
because the surface is a free surface and d{d^)/dt^V^Sd^a-^*, as 
can easily be proved by considerations similar to those in § 83. 
Thus 

47ri; =// (- uaSrd^t + VdtaSaVn - iiW ^trSd^a ;) (29). 

In the problem we have to discuss m = 0, so that iv gives only 
terms which we discussed in considering Sir Wm. Thomson's 
theory. Observing that if in w we change md% into — Sd^a we 
get for the part of w containing ?/i, w' ; we see that w' only gives 
terms that we have virtually discussed under Hicks's theory. We 
have however entirely neglected v. Are we justified in this ? In 
the first place we have seen that if the bubbles have not always 
been associated with those parts of the fluid with which they now 
are there is round every bubble absolutely no circulation. This 



* 



[Note added, 1892. This should be d (dl)ldt= -vid7:-\-^iSdZ<T^, and there 
fore equation (29) should be 

47ru = jj{ VdlaSa^u - u {(rSTdZ - inVd^<T+ I'ViO-SV/ilo-i) \ . 

This does not affect our present problem because in = in our case.] 

M. 8 



106 THE VORTEX-ATOM THEORY. [§ 98. 

shows that for any one bubble jjVd%<r = 0*, and therefore we are 
justified in neglecting the last term in equation (28). We are 
probably also justified in neglecting the first term, for probably t 
is very nearly tangential to the surface and therefore Srd^ = 0. 

[Note added, 1892. The whole of this last section is in rather 
a nebulous stage, and since writing it I have not had sufficient 
leisure to return to the matter. I hesitated whether to include it 
in the present issue. But since, notwithstanding the absence of 
any reliable results, it serves very well to illustrate how investiga- 
tions are conducted by Quaternions, I have thought it worth pub- 
lishing. 

Should anybody feel inclined to attempt to apply the method 
or an analogous one it is well to note that in the Phil, Mag. June, 
1892, p. 490, I have given the more general result sought in this 
last section. As the result is not there proved I give one proof 
which seems instructive. It exhibits the great variety of suitable 
quaternion methods of dealing with physical questions. It fur- 
nishes incidentally a fourth quaternion proof of the properties of 
vortices. It also illustrates how special quaternion methods 
developed for use in one branch of Physics at once prove them- 
selves useful in other branches. 

Adopting the notation and terminology of PhiL Trans, 1892, 
p. 686, ^ 5 — 7, let a and t be taken as an intensity and flux respec- 
tively, T still being =VV(t and therefore t = FVV. Thus er' 
and T are the actual velocity and double spin respectively and 
the equation of motion is 

Putting a = x'^(T and operating on the equation by ;j^' 

* One way out of many of proving this is as follows. Divide the bubble up into 
a number of infinitely near sections by planes perpendicular to the unit vector o. 
For any one section ^SdfX7 — 0. Consider dZ to be the element of the surface cut 
off by the following four planes, (1) the plane of section considered, (2) the con- 
secutive plane of section, (3) the two planes perpendicular to d/a, and through the 
extremities of the element dp. Thus if x be the distance between the two sections, 

dp=x~^VadI,y whence 

JJ-SadS(r = 0, 

for the surface of the bubble between the two sections. But adding, we may sup- 
pose this integral to extend over the whole bubble. Thus *SaJJ rrfS(r=0 for the 
whole bubble ; therefore a being a quite arbitrary unit vector we have for the whole 
bubble JjFrf2(r=0. 
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Now 

Xd (x"')ldt . <7 = - ;^ V-V = - X V = V,Sa,W = V . a^y2. 

.-. o- + V.o-72 = -V(i; + P) (A), 

which can easily be deduced from or utilised to prove Lord 
Kelvin's theorem concerning '* flow," equation (23) § 72 above. 

As d/dt is commutative with V, f = or t is an absolute con- 
stant for each element of matter. This being interpreted at once 
gives the well-known properties of vortices in their usual form. 

If in the equation 4i'7r(v-^P) = S.V^fJfu(v-\-P)d% we carry 
one V across the integral sign, get rid of its differentiations which 
affect u by equation (9) § 6 above, and then do the same with the 
other V we get 

47r(v + P)=fJ{(v-hP)Sd^Vu-uSdtV{v + P)}-\-fJJttV^(v+P)d%. 

At surfaces of discontinuity in o-, v and P will both be continuous, 
so that instead of f J (v -\-P) Sd^Vu we may write fjb(v +P) SdX^u. 
In the last equation substitute throughout for V(v + P) from 
equation (A). Thus 

47r {v + P)= JJb (v + P) SdXVu + ffuSdl^a - JJJuS^ad% 

+ (JJuSdXV . (7-/2 -///t.VV-ds/2) 
=JJb {v + P) SdX^ic +JJuSd%& -JJJuSV&dS -\-JJJSVuV (o-'^) ds/2, 

which is equation (36) of the Phil. Mag. paper. 

If the standard position and present position of matter coincide 
it is quite easy to prove that 

dSdX'a/dt = Sdl (& + VV^aa^) 

d (svvd%ydt =sv{&-\- vv^aa;) d%. 

Substituting for Sd^&, SV& from these in the last equation we 
get equation (32) of the Phil. Mag. paper ; but this equation can 
also be proved directly. It should be noticed that equations (34) 
and (35) of that paper have been wrongly written down from 
equation (32). In each read + V (0-72) for ~ V (0-72).] 
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